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Introduction 


I offer this work as the mathematical principles of philosophy, for the whole 
burden of philosophy seems to consist in this — from the phenomena of 
motions to investigate the forces of nature, and then from these forces to 
demonstrate the other phenomena; and to this end the general propositions in 
the first and second Books are directed. In the third Book I... derive from 
celestial phenomena the forces of gravity with which bodies tend to the Sun 
and the several planets. Then from these forces, by other propositions which 
are also mathematical, I deduce the motions of the planets, the comets, the 
Moon, and the sea, 


Sir Isaac Newton, from the Preface to the First Edition of Principia. 


The third book of Sir Isaac Newton’s monumental work Principia, entitled The System 
of the World, was published in 1687, some twenty years after he had discovered many of 
its results. In it he establishes his universal law of gravitation and uses his three laws of 
motion to predict the character of the orbits of planets and their moons. Tradition has 
it that Halley, the astronomer, Hooke, the physicist, and Wren, the architect, had 
discussed at various times in 1684 the problem of determining the orbit of a planet 
around the Sun if the force between them varied as the inverse square of the distance. 
Unable to succeed with the problem, they took it to Newton, only to find that he had 
already solved it, together with other problems of planetary motion. Encouraged 
constantly by Halley, Newton finally published Principia, The theory which it 
expounded was extraordinarily successful. Only in the twentieth century, with the 
theories of relativity and quantum mechanics, has it become evident that there are 
natural processes not explained accurately by the Newtonian world-view, 


The main aim of this unit is to apply Newtonian mechanics in order to derive Kepler's 
three laws of planetary orbits. Johann Kepler (1571-1630) arrived at these laws by trial 
and error, from careful analysis of the measurements of planetary orbits taken by 
himself and by Tycho Brahe (1546-1601), the Danish astronomer. The first two laws 
were published in 1609 and the third in 1619, so they precede Newton's mathematical 
derivations by a good lifetime, 


In deriving the laws of planetary motion we shall use modern notation and methods, 
although it must be said that the calculus in its original form was the innovation of 
Newton and of his contemporary, Leibnitz. The modern use of geometric vectors was 
apparently an invention of J. Willard Gibbs (1839-1903), and is a major notational 
advance over the geometric methods which held sway in Newton's era. 


In summary, then, there is little new mathematics introduced in this unit. Instead we 
utilize the methods of calculus and vector analysis developed in the course to describe 
one of the great scientific advances of the last three centuries, 


Study guide 


The material in Section 1 is largely revision, involving methods which will be of 
particular use later in the unit. A review of energy conservation for one-dimensional 
motion (from Unit 7) leads on to the corresponding topic for a particle moving in three 
dimensions. Then the planar motion of a particle is described in terms of polar 
coordinates and the associated unit vectors. This is a generalization of the treatment in 
Unit 28, where the particle motion was confined to a circle. 


The time spent on Section 1 will depend upon your confidence with the previous 
mechanics in the course, but you should not spend more than 30% of your total study 
time here. The major part of your effort should be reserved for Sections 2 and 3. 


Kepler's laws for planetary orbits are stated and explained in Section 2. Following this, 
Newton’s universal law of gravitation is introduced. In Section 3 the equation of 
gravitational orbits is derived, and it is shown how Kepler’s laws follow as a 
consequence of Newtonian mechanics. 
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The television programme in Section 4 covers a mathematical development similar to 
that in Section 3, but with greater reference to non-gravitational examples. More 
benefit will probably be obtained from the programme if you have managed to complete 
most of Section 3 first, but it should still be of help if viewed at an earlier stage. 


1 Review of energy and planar motion 


The ultimate goal of this unit is to model the motions of planets, satellites and comets. 
To facilitate this study it is useful to review some aspects of Newtonian mechanics, and 
that is the purpose of Section 1. 


Subsection 1.1 will remind you about the material in Unit 7 concerning particle 
motions in one dimension for which the total mechanical energy is constant. In 
Subsection 1.2 the corresponding motions in three dimensions are considered, and it is 
shown that the total mechanical energy is conserved provided that the force acting on 
the particle is conservative in the sense of Unit 26. This conclusion applies also in one 
or two dimensions, 


Subsection 1.3 shows how motion in the plane may be described in terms of polar 
coordinates r,@ and the associated unit vectors e,,e9. You saw circular motion 
expressed in this way in Unit 28. Following the derivation of kinematical formulas, 
there are some applications of Newton's second law to find equations of motion in this 
coordinate system. 


1.1 Mechanical energy in one dimension 
Consider a particle of mass m, moving in one dimension under the action of a force, We 
choose the x-axis to coincide with the line along which the particle moves. Then, as 
you know from Unit 7, Newton’s second law is mi = F, where F is the 2-component of 
the force. In general, the force may depend upon the location x of the particle, its 
velocity +, the time ¢, or any combination of these variables. We consider here the case 
where F depends only upon the particle's position «. (A special case of this is when F’ 
is constant.) A position-dependent force in one dimension is called conservative 
because, as you may recall from Unit 7 and as we shall shortly verify anew, the total 
mechanical energy is constant (conserved) for such a force. 
In terms of the force F = F(x), the potential energy U = U(x) of the particle is defined 
by the equation 
dU 
F(x) = -—, 1 
deer (1) 


or equivalently by 
U(z) = -[ F(s) ds, 
9 


where ao is any conveniently chosen point. 


Exercise 1 


(i) Find the potential energy U(r) corresponding to the force F(x) = k, where k is a 
constant, 

(ii) Find the potential energy U(.r) corresponding to the force F(x) = —kw, where k is a 
constant, 

(iii) Find the force F(x) corresponding to the potential energy U(r) = k/s (« > 0), where k is 
@ constant, 


[Solution on page 47) 


When it is possible to define a potential energy function then the sum of the kinetic 
energy 4m? and the potential energy U(r) is constant throughout the motion of the 
particle. The sum 


E=}mi? + U(x) (2) 


Unit 7 Section 3 


Unit 26 Subsection 4.3 


Unit 28 Subsection 4,2 


You may see a particle 
described elsewhere as either 
a ‘point mass’ or a ‘mass 
point’. 


The description 
‘conservative’ was not used 
when energy was introduced 
in Unit 7 Section 3. You will 
see in the next subsection 
that energy-conserving forces 
are also conservative in the 
sense of Unit 26, 


The point xo is the datum of 
the potential energy function 
U, for which U(«o) = 0. This 
point (or equivalently, the 
constant of integration 
arising from the indefinite 
integral for U) is usually 
chosen so as to achieve the 
simplest possible form for 
U(e). 
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is called the total mechanical energy of the particle, and its constancy is expressed by 
the equation E = 0. To establish that F is constant, we differentiate Equation (2) with 
respect to time t. Using the chain rule on the right-hand side, this leads to 


: ae a s » , a 
sanee Sen slnee 2), 


But the term in brackets vanishes since, by Newton's second law and Equation (1), 


Hence we have E = 0, The fact that E is constant for a conservative force often 
provides useful information about the particle motion. For example, consider the 
harmonic oscillator (that is, a particle attached to one end of a perfect spring whose 
other end is fixed), Measuring « from the equilibrium position of the particle, the force 
F is given by 


F(x) = —kz, 
where k is the stiffness of the spring. The potential energy function is therefore 
U(x) = tha? +C, 


where C is a constant. Any value of C may be chosen here since the value of 

F(x) = —dU/da is unaffected by the choice. One often takes C = 0, Figures 1(a) 
and 1(b) below show the potential energy U(s) for a negative constant C’ and for the 
more conventional choice C = 0 respectively. 


U Ui 


at 
a? 


(a) (b) 
Figure 1 


Bither form for U(x) gives the correct equation of motion. For any chosen value of C 
the total mechanical energy E = dma? + U(x) is a constant, although its value 
depends on the choice of C. 


‘To see how useful this knowledge is, suppose that initially, at ¢ = 0, the particle has 
position ro and is moving with velocity vp. This information establishes the total 
mechanical energy of the particle for all later times ¢, since E(t) = E(0). Whatever 
value of C is chosen, we have 


E= }mvj + dhaj +C = bmi*(t) + dha) +. 


Solving for the speed squared gives 
; 2 2 
a= alee pha?(t) -—C) = mE -Ula)). 


Since the left-hand side can never be negative, we know that the term in brackets on 
the right-hand side is also non-negative. Hence, for all times t we know that 


U(x) < EB. 


This information is conveyed graphically in Figure 2, where we have chosen the value of 
C used in Figure 1(a). The constant value of E is represented by the horizontal dashed 
line, At each value of x the vertical distance between this horizontal line and the curve 
for U(x) gives the kinetic energy }/nié? of the particle at that value of x. This 
difference, £ — U(«), does not depend on C. When the particle is located at x = a or 
\t x = (J, its kinetic energy is zero. The particle approaches these points, comes 
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momentarily to a halt, and then reverses its direction of motion. For this reason, points 
like a and are called turning points. 


Note that the motion of the particle is periodic: at any point x» between the turning 
points the speed is \/2(# — U(p))/m, and this value is always reproduced when the 
particle returns to that point. The velocity at the point takes the values 

+)/2(E — U(x»))/m alternately, depending on whether the particle is passing through 
®» from left to right or from right to left. The particle passes through «x, twice in each 
complete cycle. 


In fact, you know from the explicit solution for simple harmonic motion in Unit 7 that 
in this case the particle oscillates sinusoidally between the turning points. 


Exercise 2 


‘The position function for a particle of mass m, undergoing simple harmonic motion due to a 
perfect spring of stiffness k, can be written as 


a(t) = ao coswt +“ sinwt, 
w 


where w = \/k/m is the angular frequency and ro, vo are respectively the position and velocity 
of the particle at time ¢ = 0. Verify explicitly that the total mechanical energy 
E = }mi? + the? +C 
is constant, and find the turning points for this particular motion. 
[Solution on page 47} 


Explicit formula solutions to Newton's second law for the case of conservative forces 
cannot always be derived. If the potential energy function U/() is known, however, 
then the approach used prior to Exercise 2 may be applied. This gives qualitative 
information about the particle motion even when the position function can be 
established only by numerical computation. In all cases it is useful to know that EB is 
constant. Consider, for example, the motion of a particle under the influence of the 
potential energy function shown in Figure 3. 


Figure 3 
Given only the information provided by this graph, we can make the following 
statements about the particle’s motion if it has the energy F indicated in the figure. 


(i) The motion is periodic, and is confined between the turning points a and <. At 
these points the particle slows to a stop and reverses its direction of motion. 


(ii) At any point 2 for which a < x < «, the particle’s velocity is +,/2(E — U(a))/m. 


(iii) The particle is moving fastest at « = 3, since E — U(x) has an overall maximum 
there. 


(iv) As the particle passes through « = ¥ in either direction, it slows down and then 
speeds up again. This follows from the fact that the speed \/2(B — U(a))/m has a 


local minimum there. 


(v) As the particle passes through x = 6 it experiences a local maximum of speed. 


Unit 7 Subsection 2.3 
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Example 1 


What can be said about the motion of a particle under the influence of the potential 
energy function shown in Figure 4, if the energy FE of the particle is as indicated on the The potential energy 
figure? In particular, how many turning points are there at this energy level, what are function of Figure 4 is the 


these turning points and where does the particle move fastest? same as that in Figure 3, but 
the total mechanical energies 


of the particles are different. 
vu 


at 


Figure 4 


Solution 


The particle is confined (see Figure 5) to periodic motion either between points a and 
7 or between points 6 and ¢, depending upon its initial position (0), 


Figure 5 


For instance, if it were released from rest at point a then it would pursue a periodic 
motion between the turning points a and +, and its maximum speed would occur at 

x = 3. Alternatively, it might be confined by its initial conditions to move between 5 
and ¢, with maximum speed at x =. Motion at this energy E is not possible between 
points 7 and 6, as the ‘kinetic energy’ E — U(x) would be negative in this region. For 
the same reason, motion at the given energy is impossible for x <a or forx>¢. O 


It is also possible with such energy diagrams to consider the motion of a particle which 
is not bound to periodic motion. Such unbound motion occurs for the case depicted in 
Figure 6, where U(x) = 0 for x > 6. 


For example, suppose that the particle were released from rest, at 2 =a, which fixes its 
total mechanical energy for all subsequent times at E = }mz*(0) + U(a) = U(a). The 
particle would then accelerate to the right until it picked up maximum speed at x = 3. 
From z = § to « = ¥ it would decelerate to a local speed minimum at x =. It would 
then accelerate again to the right until it reached x = 6, after which it would continue 
to the right with constant speed, in unbound motion. The speed is constant for « > 6 

because the potential energy in that region is constant (zero, in fact). 
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Figure 6 


Exercise 3 


For the particle of Figure 6, released from rest at «¢ = a, give in terms of m and U expressions 
for the particle's speed at =, « = 7 and for x > 6. Arrange the speeds at the points 3, 7 
and 6 in order of magnitude. 


Exercise 4 


A particle is projected to the left in Figure 6 from some initial point zp to the right of 6, and 
with the total mechanical energy / shown in the figure, Give its initial velocity #(0) in terms 
of E, and describe its subsequent motion. 


[Solutions on page 47) 


1.2 Mechanical energy in three dimensions 


Just as for one-dimensional motion, the force acting on a particle moving in two or The gradient of a scalar field 
three dimensions can depend upon position r, velocity 7 and time ¢, or on any , denoted by grad ¢ or vo, 
combination of these variables. It often happens, however, that the force F depends was introduced in Unit 26 


only upon position, which may be denoted by writing F = F(r) or F = F(xr,y,2). If in Section, 2: It is defined as the 


addition F can be expressed as the gradient of a scalar field, then it is called a 06, . 06, . ab 
conservative force. By convention the scalar field function here is written as —U, Be! or by a ae 
where U = U(r) = U(x, y, 2) is called the potential energy of the particle. Thus in Conservative vector fields 
three (or fewer) dimensions, a conservative force F and potential energy U are were defined in terms of line 


integrals in Unit 26 


Commented yt Hel equStion: Section 3. In Subsection 4.3 


F =—gradU. (3) of that unit it was shown 
that the definition is 
Note that in one dimension this reduces to equivalent to the condition 
dU given here, involving the 
eo ab eatin gradient of a scalar field. 
ae. 


As pointed out in Unit 26 
Subsection 4,3, the potential 
energy U can be expressed as 
the line integral 


which is Equation (1) of the previous subsection. 


Conservative forces were discussed briefly in Unit 26 Subsection 4.3. It is enough for 
the purposes of this unit to recognize that where the force acting on the particle can be 
written as in Equation (3), the total mechanical energy E is constant, where - if F-dr. 


ye 
Here we have written t as a shorthand for r+, which is also equal to |#|?. In 
Equation (4) the kinetic energy bmi, the potential energy U and the total energy E 


are scalars, In order to show that F is constant, it suffices to establish that B = 0. as in Subsection 1.1. 
Differentiation of Equation (4) with respect to time t¢ and use of the chain rule produces 


(4) For the one-dimensional case 
this reduces to 


u=- f Pas, 


B= hmi* +U (x.y, 


. d 
B= (Amé-#+U(@.y,2)) ‘The derivative of a scalar 
dt product was derived in 
Sara ee olay oUiae Unit 14 Subsection 4.3. 
*  , Ov dt © dy dt ~ dz de This form of the chain rule 
=mr-i+r-gradU was introduced in Unit 25 


=f- (mi + grad). Subsection 2.4. 
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But by Equation (3), we have grad U = —F. Also F = mi by Newton’s second law, so 
the term in brackets vanishes and hence E = 0. We have shown that, when the force has 
the form given by Equation (3), the total mechanical energy is constant. This outcome 
represents the natural extension to three dimensions of the results in Subsection 1.1. 


Exercise 5 


A particle of mass m moves under the influence of a uniform gravitational field, so that 

F = —mgj, where the unit vector j is directed vertically upwards. Show that this field is 
conservative by finding the corresponding potential energy function U, 

Exercise 6 

Given that F is a conservative force with potential energy U = k/r, where k is a constant and 
r= (a? + y? + 2")"/*, show that 


5 5, kr 
(ait ui + 2k) =, 


[Solutions on page 47} 


There is a useful generalization to the result above concerning the constancy of E. 
Suppose that the particle is being acted upon by a conservative force, as defined by 
Equation (3), plus any foree G which is always directed perpendicular to the particle's 
velocity. Then its equation of motion is 


mt = —gradU +G, 


where G+ # = 0. In the exercise below you are asked to show that, under these more 
general circumstances, EF is still constant. 


Exercise 7 

Suppose that a particle is acted upon by a net force —grad U + G, where G-* = 0, Show that 
if E is given by Equation (4) then & = 0. 

[Solution on page 47] 


The result of Exercise 7 is often applied when modelling the constrained motion of 
objects. Consider, for instance, a particle of mass m sliding without friction down an 
inclined plane (see Figure 7). Since we are assuming that there is no friction present, 
the only forces acting are gravity and the reaction force R. The latter is always 
perpendicular to the plane and hence also to the velocity, The gravitational force —mgj 
is conservative (as shown in Exercise 5) with potential energy function 


U=mgy+C, 


where C is any constant, If u is the speed then, from the result of Exercise 7, the total 
mechanical energy F is constant, where 


Figure 7 


= 2 
B= xmu* + mgy + C. You saw another example of 
this type, featuring a force 
Exercise 8 directed at right angles to 
In the situation just described, the particle is released from rest at a vertical height h above _the velocity, when studying 


pee ae : ae the motion of a pendulum in 
7 ax! . ts e 
he x-axis in Figure 7. Find its speed uo when it reaches the origin Fa apeay ete er 


[Solution on page 48 | 


1.3 Planar motion using polar coordinates 


In this subsection we consider a means of describing the motion of a particle in a plane. 
For some applications it is convenient to describe the position vector r of the particle in 
terms of Cartesian coordinates (x,y) and Cartesian unit vectors i,j, that is, 


r(t) = 2(t)i + y(t)j. (5) 
Then the velocity and acceleration vectors are respectively 

v(t) = F(t) = &(t)it H(t)j (6) 
and a(t) = ¥(t) = #(t) = #(t)i+ H(t}. (7) 


The ease of derivation of Equations (6) and (7) from Equation (5) is due to the fact 
that i and j are constant vectors throughout the motion of the particle. In summary, 
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then, the motion of a particle in a plane can be described by giving its Cartesian 
coordinates (w,y) as functions of time, together with the corresponding fixed 
orthogonal unit vectors i and j. 


In Unit 28 you saw that it is sometimes more convenient to use polar coordinates [r, 4] 
and the corresponding orthogonal unit polar vectors e, and ey (see Figure 3). For 
instance, corresponding to Equation (5), the polar form of the position vector is 


rep, (8) 


where r = |r| is the distance of the particle from the origin. This equation is not as 
simple as it looks because, for general motion in the plane, both r and e, depend on 
time, as do the polar angle @ and its associated unit vector ey. 


To see this time dependence, consider the fact that any vector in the plane can be 
expressed in terms of either the pair i,j or the pair e,,e9. In particular, each of i and j 
themselves can be expressed as a linear combination of e,,es, and vice versa, From 
Unit 28, or by noting that the unit vectors e,,e9 are obtained by rotating i,j 
respectively through an angle @ (see Figure 9), these relations are 


e, = cosdi+sin0j, 
i= cose, —sin@eg, 


eo = —sin#i+ cos0j (9) 
j=sin@e, +cosep, (10) 


Equations (9) show that e, and eg are functions of the polar angle @. Thus, if a particle 
is moving along some curve so that @ is time-dependent, then e, and eg depend on time 
too, For instance, Figure 10 shows the curve followed by a particle moving anticlockwise 
around and away from the origin, which passes through the point [r,0;] at time t; and 
the point [r2,62] at some later time ft. Note that, while e, and eg depend upon 0, they 
are independent of the radial distance r (which does not appear in Equations (9)). 


and 


ve 


©, (3) 


ay 


Figure 10 


Exercise 9 


Show explicitly from Equations (9) that, at each point other than the origin, e,.e9 form an 
orthogonal pair of unit vectors, that is, they satisfy the conditions 


e--e9=0 and e,+e,=ey+ey = 1. 
(You will need to use the fact that i,j form an orthogonal pair of unit vectors.) 


[Solution on page 48) 


In order to study the motion of a particle in plane polar coordinates, it is necessary to 
have expressions for its velocity and acceleration vectors in terms of r,0,e, and eg. In 
the course of deriving these equations, we need to express the time derivatives é,,é9 in 
terms of e,.ey. The required relations for achieving this are 


@-=6es and é=—be,, (11) 


which may be obtained by differentiating Equations (9) and applying the chain rule. 


Unit 28 Subsection 4.2 


The polar coordinates [r, 4] 

are defined in terms of the 

Cartesian coordinates (:,y) 

by the equations 
a =reosé, 


y=rsind. 


Figure 8 


Equations (9) and (10) apply 
at all points of the plane 
other than the origin, where 
@,ey and 9 are undefined, 


Figure 9 


‘These equations were derived 
in Unit 28 Subsection 4.2. 
Although attention was 
restricted there to circular 
motion, the derivations apply 
also in the general case. 
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We shall now apply Equations (11) to obtain expressions for the velocity v = and 
acceleration a = 7 in terms of polar coordinates and unit vectors. As pointed out 
above, the position vector r is given by 


r=re), (8) 
where r = |r|, Differentiating this equation and using the product rule, we have 

r=fe, +réy, 
which by the first of Equations (11) may be written as 

b=fe,+rbep. (12) 
Differentiation of this expression for the velocity produces the relation 

# =e, +e, + (76+ rd)eg + rb6p 


for the acceleration. After applying each of Equations (11) and collecting terms, we 
arrive at 


# = (i —r6)e, + (270 + rd)ep. 
Recognizing that 
i cH Liens 
270 +70 = eae 0), 


we can write the expression for the acceleration as 


#=(F-ré)e, + 14 (0)e6. (13) 
The following box summarizes the results which we have obtained, 
Position, velocity and acceleration in plane polar coordinates 
‘The plane polar unit vectors e,,eg are related to the Cartesian unit vectors 
i,j by the equations 
e, =cos@i+sin#j, eg =—sindi+cosAj (9) 
and i= cose, —sinfep, j=sin@e, + cos ey. (10) 
The derivatives é,,é9 are given in terms of e,,ey by 
@, =Gey and é = —de,. (11) 
In terms of r, @, e, and eg, the position r, velocity ¢ and acceleration # of a 
particle moving in the (2, y)-plane are given by the equations 
r=re,, (8) 
t=fe, +rdey, (12) 
f= (F —16%e, 4+ 14 Heo. (13) 


Of the following four exercises, the first three will give you practice in expressing a 
particle's velocity and acceleration in plane polar coordinates. The fourth, Exercise 13, 
asks you to express the mechanical energy of a particle in polar coordinates. The result 
of this exercise will be applied later in the unit. 


Exercise 10 


A particle is constrained to move in a circle of fixed radius r = R whose centre is at the origin. 


Find the velocity # and acceleration F in terms of R, @, @, e- and eg. Give the velocity and 
acceleration also in terms of R, @, 6, 6, i and j. 


Exercise 11 


A particle is constrained to move radially outwards from the origin, at some fixed angle @ = @. 


Find the particle's velocity # and acceleration ¥ in terms of 7, # and e,. Give the velocity and 
acceleration also in terms of 7, #', 9, i and j. 


Note that # =djr|/dt, the 
rate of change of the radial 
coordinate, is not in general 
the same as the speed |r| of 
the particle. 


‘These expressions for the 
velocity and the acceleration 
differ from those derived in 
Unit 28 Subsection 4.2 
because there we restricted 
ourselves to circular motion 
with r = constant. 


This is the particular type of 
planar motion which was 
considered in Unit 28 
Section 3 and Subsection 4.2. 
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Exercise 12 

A particle moves in the (2, y)-plane with radial coordinate r(t) = at® + bt and angle 0(t) = ct”. 
Find the velocity and acceleration in terms of ¢, e, and eo. 

Exercise 13 


It was shown in Subsection 1.2 that if a particle moves under the influence of a conservative 
force F = —gradU;, then its total mechanical energy B = }mi? +U is constant. Show that if Recall that by #? we mean 
the motion takes place in the (.r,y)-plane then B may be written as Fer=|i/*. 


B= im(i? +16") +U. 


{Solutions on page 48 | 


Thus far we have derived kinematical formulas for the planar motion of a particle, 
using polar coordinates. Before proceeding to the subject of planetary orbits we shall 
briefly examine the mechanics involved in some simpler systems. The starting point for 
this examination is, as usual, Newton's second law. 


In plane Cartesian coordinates, Newton's second law mi = F takes the form 


m(#i + ij) = Fri+ Fi. Here F, and Fy are 
’ : i : respectively the a- and 
leading as in Unit 15 Section 2 to the two scalar equations y-components of the force F. 


mi =F, and mij = F,. 


Similarly, Newton’s second law mi = F may be expressed in plane polar coordinates, 
using Equation (13), as 


F ter 
m ( — re, + 2 (r*a)ee) = Fre, + Foes, Here F, and Fy are 
rdt respectively the radial and 
‘at a 5 transverse components of the 
giving the pair of scalar equations ' force F. 
m(it—ré?) =F, and (96) = Fy, (14) 


It is these equations of motion which are sought in the following example and exercises, 
each of which represents a situation in which it is appropriate to use polar rather than 
Cartesian coordinates. 


Example 2 


A particle of mass m is launched at the bottom of a circular track of radius R (see 
Figure 11). As the particle moves along the track it is subjected to the downward force 
of gravity, a normal reaction force of magnitude F perpendicular to the track, and a 
frictional force of magnitude f in the direction opposite to that of the particle’s motion. 


Set up the equations of motion for the particle in plane polar coordinates, assuming 
that the particle remains in contact with the track and moves in the sense of 
increasing 0. 


Solution 


The forces have been drawn on Figure 11, which also indicates our choice of Cartesian 
coordinate system. The reaction exerted by the track on the particle is —Fe,, The 
force of gravity is mgi, but i is given in terms of e,, and eg by the first of 

Equations (10), so the gravitational force may be written as Figure 11 


mag(cos #e,. — sin 8 eg). 


The particle moves in the sense of increasing 0, that is, in the direction of ey, so the 
frictional force is —feg. Adding the three forces together gives the total force 


— Fe, +mg(cos#e, — sin e9) — fea = (—F +mgcos6)e, + (—f — mgsin®)en, 
so that the radial and transverse components of the total force are respectively 


F,=—F +mgcos@ and Fy =—f —mgsind. 
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Since the particle remains in contact with the track, it is constrained to move in a circle 
with constant radial coordinate r = R, so that + = 7 = 0. The equations of motion 
(Equations (14) are therefore 

— mR6? = —F + mgcos0, 

mR6=—f—mgsind. 2 


Exercise 14 


A puck of mass m moves without friction around the origin O on a horizontal table (see 
Figure 12). It is fixed to a light, inextensible string which is reeled in steadily, so that its 
distance from the origin is given by 


r(t)=r(0)—pt — (t < r(0)/n), 
where 4 is a positive constant. The string remains taut throughout. 
(i) Show that the equations of motion for the particle which represents the puck are 
di; a, 
a (778) =0, 
where T is the tension in the string. 


mré?=T and 


(ii) Show that at time ¢ the tension is given by 
mr*(0)62(0) __ mr*(0)62(0) 
rat) (r(0) = pt)” 
(iii) Suppose that m= 0.1kg, 1 = gj ms”, the initial distance is r(0) = 1m and the initial 


angular velocity is 6(0) = 1rads~', Suppose further that the string will break when it 
experiences a tension of 500N, Find the distance 7 and the time ¢ at which the string 
breaks. 


T(t) = 


Exercise 15 


A puck of mass m moves frictionlessly around the origin on a horizontal table. It is connected 
to the origin by a perfect spring of natural length lo and stiffness k, so that it is acted upon by 
a force 


F= 


(i) Set up the equations of motion for the particle which represents the puck, 


(r = lo)er. 


(ii) Show that, however the particle moves, mr?4 is constant, 
(iii) 


(iv) 


Verify that a potential energy function for this force is U = 4k(r —Io)?. 


Suppose that the puck has speed u when the spring has length 2/9. What is its speed 
when the spring has its natural length in the subsequent motion? 


[Solutions on page 48 | 


Summary of Section 1 


1. For particle motion in three or fewer dimensions, a force F acting on the particle is 
conservative if there exists a scalar field function U = U(r) = U(a,y, z) such that 


F =~gradU. 


The function U is the potential energy of the particle, In such a case, the total 
mechanical energy of the particle, 
E= hme? +U 
(where # = ¢- # = |#|?), is constant throughout the motion. 
2. In one dimension, the condition for a conservative force and the expression for the 
total mechanical energy become respectively 
du 


di: 
where F is the z-component of the force and U = U(x). 


and E = imi? +U, 


With f = 0, these are also. 
the equations of motion for a 
circular pendulum, where F 
is now the tension in the 
string. You saw these 
equations derived via 
Cartesian coordinates in 
Unit 15 Section 5, and more 
directly in Unit 28 
Subsection 3.3. 


‘The upper bound on the time 
t is included because 7(t) 
must always be non-negative. 
The puck would reach the 
origin at time t = r(0)/j if 
the situation remained as 
described until that time. 


Figure 12 


You will see this apparatus 
demonstrated and discussed 
further in the television 
programme for this unit 
(Section 4). 
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3. The unit vectors e,,eg corresponding to the plane polar coordinates [r, 4] (where 
a“ =rcos@, y =rsin@) are related to the Cartesian unit vectors i,j by the 
equations 

e, = cos@i+sin@j, eg = —sin@i+cos0j 
cose, —sin@ eg, j=sin@e, +cos0e,. 


and 
4. The derivatives é,,é are given in terms of e,,eg by 
@,=0eg and é = —ée,. 


5. The position r, velocity ¢ and acceleration # of a particle moving in the 
(a. y)-plane are given in terms of plane polar coordinates by the equations 


r= Tey, 
f= fe,+r0eg, 


Sie My ce 
F=(7 —70' Jen + = alr" Aes. 


2 Kepler’s laws and Newton’s law of 


gravitation 
The Earth is an ordinary-sized planet which orbits a middle-aged star of rather small 
size near the outer fringes of a typical spiral galaxy some 10° light-years across. This One light-year is the distance 
galaxy contains between 10'' and 10!* visible stars, together with much other matter, travelled by light in one year, 
and is one member of a rough association, or cluster, of some thousand galaxies of or approximately 
various shapes and sizes. a 
Located in this myriad of astronomical objects is our solar system. It consists of a 
relatively massive Sun encircled by eight chief planets, some of which themselves Owing to its small size and 
possess one or more satellite moons. Any member of the solar system interacts via features of its orbit, Pluto is 
gravitational force with all of the others and, indeed, with all other heavenly bodies. ears iconsidered: to 


But distances between stars are measured in light-years—and the force of gravity 
between two bodies decreases rapidly in magnitude as their separation increases—so 
interactions with objects beyond the solar system can mostly be neglected. 


The gravitational influence of a body increases with mass, and within the solar system, 
there is a rough hierarchy of size. The Sun is by far the largest object, followed by the 
planets and then by their moons. The Earth’s mass is only about 3 x 107° that of the 
Sun, and its moon's mass is about 10? of its own. The largest planet, Jupiter, is 
roughly 300 times more massive than the Earth, but it is far away. In fact, the 
gravitational influence of a planet upon any other body is of little significance 
compared with the effect of the Sun, unless the body acted upon is relatively close to 
the planet. To a good first approximation, then, we can think of the planets as orbiting 
the Sun singly, and the moons as singly orbiting their planets, This is the model 
adopted by Newton, who considered the problem of only two bodies orbiting one 
another. (In a more advanced treatment, the effects of the other planets can be added 
as perturbations.) We shall assume in addition that one of the bodies is very much 
more massive than the other. 


An analysis of the motion of one relatively light body moving under the gravitational 
attraction of a second, relatively massive, body is the goal of Section 3. In this section 
we introduce the topic of planetary motion by describing (in Subsection 2.1) the 
features of their orbits which were observed by Kepler and stated in the form of three 
informative laws. These laws will be derived in Section 3 as a consequence of Newton’s 
universal law of gravitation, which we consider in Subsection 2.2. 
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2.1 Kepler’s laws described 


As they are worded, Kepler's three laws refer specifically to the orbits of planets around 
the Sun, but they apply also to the orbit of any relatively light body about another 
much more massive body, The laws were established by Kepler from experimental 
observations of the motion of several of the planets. You will see later that they may be 
derived from Newton's three laws together with his law of gravitation. It is the purpose 
of the current subsection to explain the meaning of Kepler's laws. 


Kepler’s laws of planetary motion 
Law I Pach planet moves in an ellipse, with the Sun at one focus. 


Law II The line joining a planet to the Sun sweeps out equal areas in equal 
times. 


Law III The square of the orbital period of a planet is proportional to the 
cube of the semi-major axis of its orbit. 


Kepler’s first law specifies the shape of a planetary orbit. Since it refers to an ellipse 
and to a focus, we start by considering the mathematics of ellipses. 


An ellipse is a symmetrical plane figure like that drawn in Figure 1. One of the 
implications of this first law, then, is that the orbit of a planet about the Sun lies in a 
plane. The usual definition of an ellipse is as a set of points (X,Y) which, with a 
suitable choice of origin and axes, satisfy the equation 

x? y? 
ptpah (1) 
where a and 6 are positive constants with a > b. The constants a and b are called 
respectively the semi-major ais and semi-minor aris of the ellipse. Their geometrical 
significance is shown in Figure 1. 


Exercise 1 


What can be said about an ellipse in the special case a = b? 


Exercise 2 


(i) Show that the area A of an ellipse can be expressed as 


A= a f° 1— (X/a)* dx. 


(ii) By using the substitution X =a cos, show that A = mab, 
[Solutions on page 49 | 


To understand Kepler's first law, it remains to explain where a focus of an ellipse is 
located, and then to express the equation of the ellipse relative to this point. We shall 
express the equation with respect to the new set of axes Oxy shown in Figure 2 below. 
These axes are obtained from the set O'XY of Figure 1 by a translation of the axes to 
the right by a distance g, where 


q=Va—8. (2) 


(We could also shift axes to the left by the distance g, but our choice represents no loss 
of generality since this configuration is just the other rotated by an angle 7.) The 
points F and O of Figure 2, located a distance q to the left and to the right of the 
centre of the ellipse, are called the foci of the ellipse. 


Figure 1 


‘Foci’ is the plural of ‘focus’. 
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Figure 2 


Exercise 3 
Where are the foci of a circle? 


(Solution on page 49 | 
Since the translation is along the coincident «- and X-axes, the new coordinates (x,y) 
are related to the original ones (X,Y) by 
X=2+q, Y=y 
‘Then, from Equation (1), the equation of the ellipse takes the form 


e+q)y . y? 

(x =a) ah 

a bo 

We shall now proceed to re-express this equation in plane polar coordinates. After 
multiplying the equation by 6? and rearranging, we obtain 


=1. 


2p? 2p Pes wera) 
yao — Beta ee — aie" age he): 
Adding xx? to both sides and collecting terms gives 
v? 2b7q a 
Caray rent (jt Bema 2) Va see a, 
oy =(° at) we (1 a). 


The left-hand side is r?, the square of the radial distance, so that, on using 
Equation (2) to express g in terms of a,b, we have 


ny 2 2 2 
ieee Pa 1-3+(1-5) 2. (3) 
a 


a a? 


This rather cumbersome expression is made neater by employing two new quantities 

defined in terms of a and b. These quantities are the eccentricity This usage of the symbol e 
should not be confused with 
its use as the base of natural 


4 
(4) jogarithms, 2.718 28. 


@=-= 


a 
and the semi-latus rectum 
e 
=—. 5 
= (5) 
Note that, since b < a, we have 0 < e < 1. This is a condition on the eccentricity for 
ellipses. 
Exercise 4 
Find e when 
(i)a=6 (a circle); (ii) b= $a; (iii) b= fa. 


Comment on the relationship between the eccentricity and shape of an ellipse. 


[Solution on page 49) 
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From Equations (3)-(5), we now have 
r? =P — lew + ea? 
or r? = (l—ex)*, 
Then taking the square root of both sides gives 
r=+(l—ez). (6) 


We must now decide on the choice of sign for the right-hand side of Equation (6), to be 
consistent with the fact that r is never negative. We shall show that, for all points x on 
the ellipse, the quantity / — ea is positive. From Figure 2, the largest value of x on the 
ellipse is 7 =a —q. Thus we have, for all points on the curve, 


l—ex >l—e(a-q). 
But we also have q = Va? — 6, | = b?/a and e = q/a, so the right-hand side of this 
inequality is 


2 Vara 
1-e(a—q)= £ S we (a- v=) 


2 a_i 
=~ van—e ye —4 
=a-Va—B, 


which is never negative. Hence | — ex: itself is never negative, when 2c assumes values for 
points on the curve. Thus the positive sign is appropriate in Equation (6), and the 
equation for an ellipse with respect to the axes Owry of Figure 2 is 


Since x = rcos@ (see Figure 3) we have the convenient form 
l 
= 1+ecosé (7) 


for the equation of an ellipse in polar coordinates with origin at a focus of the ellipse. 


Kepler's first law therefore says that if the origin of coordinates is located at the Sun, 
then the orbit of each planet lies in a plane containing that origin, and is described by 
Equation (7) once the axes have been appropriately oriented. With a different 
orientation for the axes (but the same origin), Equation (7) becomes 


! = 1+ ecos(4 — 0), (8) 


where 9 is the clockwise angle through which the axes have been rotated from their 
position in Figure 3. The major axis of the ellipse then lies along the direction specified 
by 0 = %. 


Kepler’s second law says that the line joining a planet to the Sun sweeps out equal 
areas in equal times. This law can be used to specify the rate at which a planet moves 
in its particular elliptical orbit. 


Suppose, without loss of generality, that the planet moves in an anticlockwise sense. In 
Figure 4 the two shaded regions represent the areas swept out by the planet between 
times ¢, and ta, and between times ts and ty. Assuming that the two time intervals are 
of equal duration, Kepler's second law says that these areas are equal. If A(t) is the 
total area swept out since some initial time, ¢ = 0 say, this law can be expressed as 


A(ta) — A(ts) = A(t) — A(ti) whenever ty — ts = ta — th. 


Putting « = 0 into this 
equation gives r = +1, 
showing that the semi-latus 
rectum / is the distance from 
the origin of the ellipse’s 
intercepts on the y-axis. This 
is illustrated in Figure 3 
below. 


Figure 3 
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at 


Figure 4 


A more concise way of stating Kepler's second law is that the rate at which area is 
swept out is constant, that is, 


A is constant. 


It is possible also to state the law in a form which relates the radial coordinate r to the 
rate of change @ of the angular coordinate. Figure 5 shows a close-up of part of the 
planetary orbit, where we consider the area 6A swept out during a small time interval Eo 
of duration 6t, Over this interval the planet moves from point P at time t to point P’ 

at time t + 6t, traversing an angle 60 and changing its radius from r tor +ér. Ifétis Figure 5 
sufficiently small then the area 6A will lie between the areas of two circular sectors with 


angle 60, one having radius r and the other radius r + 6r. If r(t) is increasing in the If r(t) is decreasing then the 
part of the curve under consideration then this gives expressions for the upper and 
j : lower bounds in this 
$7760 < 6A < 4(r + 6r)*60. inequality are reversed, 


A fed P leading to the same outcome. 
On dividing through by 6t, and then taking the limit as 6¢ tends to zero, we obtain 


A= 1r76. (9) 
Since A is constant according to Kepler's second law, this law may also be interpreted 
as saying that 

r76 is constant. 


Equation (9) can also be used to express an area swept out by the planet in terms of r 
and @. Thus the shaded area in Figure 4 covered between times t; and ty is 


ta 
Atta) ~ ath) = [ale = f Aat= [ 1,2 (1)6(0) dt. 


Exercise 5 


A body moves in a circular orbit around the origin, with radius R and fixed speed V. Show 
that A is constant, and find its value. 


Exercise 6 


Compare the average speeds of the orbiting planet during the two equal intervals t; <t < ta 
and ty <¢ < ts in Figure 4. Your answer should consist of a few sentences only. 


[Solutions on page 49] 


Kepler’s third law gives a relation between the period of the planetary orbit and the 
lateral extent (semi-major axis) of the orbit. If T is the period then, by definition, the 
planet will take this time to make one full circuit of the ellipse. If the particular orbit 
has semi-major axis a, then Kepler’s third law says that T? is proportional to a%, 

that is, 


Tra", 
where / is a constant. This law, then, makes a statement connecting the overall rate of 


motion of the planet to a geometrical property of its orbit. We shall prove this result, 
together with Kepler's other two laws, in Section 3. 
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Exercise 7 


The semi-major axis of the planet Mars is 1.524 times as large as the Earth’s semi-major axis. 
The Earth's period for its orbit around the Sun is 365.256 days. Find the period (in Earth 
days) for the orbit of Mars around the Sun. 


[Solution on page 49] 


2.2 Gravitational force and potential energy 


I began to think of gravity extending to ye orb of the Moon, and ... from 
Kepler’s Rule ... I deduced that, the forces which keep the Planets in their 
Orbs must [vary] reciprocally as the square of their distances from the centres 
about which they revolve: and thereby compared the force requisite to keep 
the Moon in her Orb with the force of gravity at the surface of the Barth, and 
found them to answer pretty nearly. All this was in the two plague years of 
1665 and 1666, for in those days I was in the prime of my age for invention, 
and minded Mathematicks and Philosophy more than at any time since, 


Isaac Newton's account of his discovery of the universal law of gravitation. 


This subsection starts with the statement of Newton’s universal law of gravitation, then You met Newton's law of 
derives the potential energy function corresponding to the (conservative) gravitational gravitation in the television 
force exerted by a particular type of extended object. Finally, we compare the theory of Programme for Unit 28, 
gravitation introduced here with the simplified version which has been applied in 

previous units for motion near the surface of the Earth. 


Newton’s universal law of gravitation describes the nature of the force between pairs of 
particles. It gives the attractive force between them as a function of the masses of the 
particles and of their separation. It is universal in the sense that it is hypothesized to 
hold between the particles comprising all bodies in the universe, and not just between 
those comprising the Earth and an apple, or the Earth and the Moon, In words, the 
law can be stated as follows. 


Newton’s universal law of gravitation 


The force of gravity between two particles is attractive, directed along the line 
between them, proportional to the product of their masses, and inversely 
proportional to the square of their separation. 


Let us write this mathematically, using the vector notation which was unavailable to 
Newton. We write the position vectors of the two particles as rj,rj and their masses as 


mj,mj;. The gravitational force on m, due to the presence of m j is denoted by Fy;. As a shorthand device, we 
From the verbal statement of Newton's law above, we know the following. shall refer to the particles by 
using their masses mj, my as 
(i) Fi; is directed from m; to mj, hence its direction is defined by the unit vector labels. 
Li 


(ii) The magnitude of F,, is proportional to mymy, and to ent 
4 Ti 


Putting these elements together gives the force of gravity on m, exerted by m; as 


= ——— (ry — 1): (10) 


This is Newton’s universal law of gravitation in vector form. Here G is the 
constant of proportionality, called the universal constant of gravitation, or gravitational Figure 6 
constant. It is one of the least precisely measured physical constants. Its value, in SI 

units, is 


G ~6.67 x 10-". (11) 
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Exercise 8 
What are the SI units of G? 


Exercise 9 

Show that © = 5 is a unit vector. 
7 

Exercise 10 


Gravitation is an inter-particle force. By Newton’s third law (Unit 17 Section 2), the force 
exerted upon mi by my should be equal in magnitude, but opposite in direction, to the force 
exerted upon m, by m;. Show that the universal law of gravitation is consistent with this, so 
that Fy, +F) =0. 


[Solutions on page 49 | 


Most of the rest of this subsection is taken up with consideration of the gravitational 
attraction caused by an extended object with spherically symmetric mass distribution, 
and calculation of the corresponding potential energy function. 


Using Equations (10) and (11), it is possible to calculate the force of attraction between 
any two objects, provided that they are sufficiently small to be modelled by particles. 
In many cases, however, it is far from the truth to claim that both objects are ‘small’ 
relative to the situation being considered. For example, Figure 7(a) shows a satellite of 
mass m., in orbit around the Earth at a height approximately equal to the Earth's own 
radius. At first sight it seems unreasonable here to model the Earth by a particle, 
although it appears sensible (for most purposes) to use a particle model for the satellite. 


am am 


(a) (h) 
Figure 7 


In fact, you will see in this subsection that if the Earth’s mass is assumed to be 
distributed in a spherically symmetric manner (to be defined below), then its attraction 
on the satellite is the same as if all its mass were concentrated in a particle at its centre 
(see Figure 7(b)). 


To approach this problem, consider the gravitational force on some particle of mass m. 
outside an arbitrary extended object. We want to calculate the net force acting on the 
particle with respect to some coordinate system Oxyz, as in Figure 8. To do this, we 
regard the extended body as being composed of many small volume elements. Then the 
net force on the exterior particle is obtained by adding together all the (vector) 
contributions from these volume elements. Since the elements are small, they can be 
treated individually as particles. 


If you are short of time, you 
might like to read as far as 
the paragraph below Figure 7 
and then move directly to the 
summary box on page 25. 
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Figure 8 


The force exerted upon the particle of mass m by a typical volume element, which has 
mass m; and position vector rj, is obtained by applying Equation (10). This force is 


r-r 
Gmm,— 


where r is the position vector of the mass m. Adding together the contributions from all 
the volume elements in the extended object gives the net gravitational force on m as 


mi(r; — 1) 


N 
F(r)=Gm yo (12) 
i= 


Ini —r/8 
By construction, the sum in Equation (12) supposes that the extended object consists 
of a large number, N, of small elements m,. We shall ultimately consider (for a 
particular case) the limit of this sum as N becomes large and each m; tends to zero, 
which leads to an expression for F(r) as an integral over the volume of the extended 
object. This could be done by taking the limit directly from Equation (12), but it is 
simpler to rewrite it first in terms of gravitational potential energy. 


You saw in Subsection 1.2 that a force F(r) is conservative, with potential energy 
function U(r), provided that. 


F = —gradU. 


We shall now show that the net force given by Equation (12) can be expressed in this 
way. The main step in this demonstration is to establish that 
yr 1 


a A —. 13 
rip j 18) 


In order to prove Equation (13), it suffices to use Cartesian coordinates, in terms of 
which we have 


a2 


Ir — ri] = ((a — ay)? + (y — ys)? + (2 — 4)?) 


The definition of the gradient gives 
i a 1 owt Od 


‘ad =f j Fi 
see =tl Oeh—n] *8yr—ed * der =n 


(14) 


The i-component on the right-hand side of this equation is 


= (2 =)? + (yu)? +@-a)) 0? 


= i(—}) ((@ — 4)? + (y—w)? + (2-2?) °? x A@—m) 
=p 2s 


- rari 


The calculations for the j- and k-components of Equation (14) are similar. Using the 
relations 


reait+yj+2k and r;=aji+ yj +k, 


‘This is the gravitational force 
on the particle of mass m. 
due to the system of NV 
particles. Although this is 
only an approximation to the 
gravitational effect of the 
extended object, the 
approximation becomes exact 
in the limit as N is increased 
indefinitely. 


Equation (3) of Section 1 
(page 9) 
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we obtain 


i} 


— (ie — mi) + i(y — 4) + k(z — 21) neaate 
Nor 


jr—ry 


which proves the validity of Equation (13). 


Combining Equations (12) and (13) then gives the gravitational force F in terms of a 
L 
Ir—ri| 


potential energy function, as 
mi 
fr ri 
my 
or F(r) = — grad (-omy> rs ) : 
iag  S 


This equation expresses the net gravitational force exerted by a system of N particles, 
with masses m, and corresponding position vectors r; (i = 1,2,...,.N), upon an 
external particle of mass m and position vector r. It has the form 


F(r) = — grad U(r), 


where the potential energy function U is given by 


N 
U(r) = -Gm>> 
i=l 


N 
F(r) = Gm }> mj grad ——— = grad (cm 
i=l 


SE (15) 


Exercise 11 


A particle of mass M is fixed at the origin. Find the corresponding gravitational potential 
energy function U for a second particle which has mass m and position vector r. Find also the 
gravitational force exerted upon the second particle by the first. 


[Solution on page 49) 


We seek now to progress from the gravitational potential energy of a particle due to a 
system of particles to the corresponding potential energy for a continuous distribution 
of mass, as in Figure 8. To achieve this, we characterize the distribution of mass by its 
mass density function, p(x). This function is defined so that any small volume element ‘The mass density function p 


having volume 6V; and centred at position r; contains a mass m,;, where was introduced in Unit 27 
Section 2, where it was used 
my ~ plr,)6Vi. to calculate the overall mass 


i : fee ‘ oaks a of extended bodies. 
Inserting this expression into Equation (15) gives, for the gravitational potential energy 


U of a particle due to a body with continuous mass distribution, 

N 

ry dV; 
Une son se eee 

a al 

ist 
This approximation to U becomes exact in the limit as the size of the volume elements 
5V; tends towards zero and their number, N, increases indefinitely. In this limit, we 
obtain 


(16) 


The right-hand side here is a volume integral with respect to the variable r’, extending 
over the entire body B. 


Equation (16) gives quite generally the gravitational potential energy of a particle of 
mass m located outside any continuous distribution of mass, specified by the mass 
density function p. We now specialize to the case of a mass m lying outside a 
spherically symmetric distribution of mass. 


By definition, a mass distribution is spherically symmetric if the corresponding mass 
density function p(r) depends only upon the radial distance, r = \/x? + y? + 22. We 
express this spherical dependence of the mass density function by writing it as 


alr) = e(|r|) = p(r). 
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Exercise 12 


Which of the following mass density functions p correspond to spherically symmetric mass 
distributions? (A, B and R are positive constants.) 


(a) p(r)=A(n?+y'+2%)e-B" — (r > 0) 
(b) p(r)= A(x? +y?)e-B" (rn > 0) 
(©) p(r)=A(e?+y?+24)?e-8" (30) 


: 
@ oat gr Psts) 


si os 
©) ote) = {ge aS reac 


© o@=49 PSRr™ 


[Solution on page 49 | 


The integral on the right-hand side of Equation (16) was considered in Unit 27 for the 
case of a spherically symmetric mass distribution. It was proved there that when 
A(t’) = p(|r'|) then 


Fi A ele (a7) 
'p |r! — | i 
where 
M= f oe)av (18) 
B 


is the total mass of the sphere B. By Equations (16), (17) and (18), the gravitational 
potential energy function for a spherically symmetric mass distribution is 
U(r) = oa (19) 


Exercise 13 
Show that the gravitational force exerted on a particle of mass rm due to a spherically 
symmetric mass distribution of mass M with its centre at the origin is 


F(x) = Gm 


where e, = r/|r| is the unit vector in the outward radial direction. 


[Solution on page 49] 


Equation (19), together with the results of Exercises 11 and 13, establishes the truth of 
the statement following Figure 7: the gravitational attraction of a body with spherically 
symmetric mass distribution on an exterior particle is the same as if all of the body's 
mass were concentrated at its centre. In proving this result we have derived expressions 
for the gravitational force and potential energy which will be applied in Section 3 to the 
task of specifying planetary orbits. These expressions are given once again below. 


Unit 27 Subsection 4.3 


If the body is of finite extent 
and has a spherically 

symmetric mass distribution, 
then it is necessarily a sphere. 
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The gravitational effect of a spherically symmetric body 


A body has a spherically symmetric mass distribution if (with the origin of 
coordinates at the centre of the body) its mass density function p(r) depends 
only on the radial distance |r| =r = \/x? + y? + z?. The gravitational force 
exerted by such a body B on an exterior particle with mass m and position 
vector r is given by 

GmM 


yr ers 


F(r) = —gradU(r) = ~ 


where e, = r/r is the unit vector in the outward radial direction, 


Ug=- Sul 


(19) 
is the gravitational potential energy of the particle, and 
M= [ p(r’) dV (18) 
JB 


is the total mass of the body. 


Exercise 14 


Two spherically symmetric bodies have masses M and m. If their centres are separated at 
some instant by a distance R, what are then the magnitude and direction of the gravitational 
forces exerted by each body upon the other? (Hint: By a result of Unit 17 Section 2, the 
motion of the centre of mass of a body is the same as that of a particle of the same total mass 
which experiences all of the external forces applied to particles of the body.] 


[Solution on page 50 


To conclude this section, we examine how Newton's universal law of gravitation matches 
up with the simpler model used previously for motion near the surface of the Earth. 
When modelling motion close to the Earth it is customary to assume, as we have done 
in earlier units, that the force on a particle of mass m is directed vertically downwards 
and has the constant magnitude mg, where g = 9.81 ms~*. More generally, under the 
assumption that the Earth is a spherically symmetric distribution of mass, we have just 
shown that at any point on or above the Earth’s surface the gravitational force is 
F(r)=- mMG 


or) 


where M is the mass of the Earth, e, = r/r is the unit vector in the outward radial 
direction and r = |r| is the distance from the centre. 


Clearly the two descriptions agree as to the direction of the gravitational force, since 
~e,, being directed towards the centre of the Earth, is always ‘vertically downwards’ 
from a local point of view. It remains then to compare the two expressions for the 
magnitude. 


Exercise 15 

(i) Show that g = MG/R®, where R is the radius of the Earth. 

(ii) Given that g = 9.81ms~*, G = 6.67 x 107! m*kg~'s~? and R = 6.4 x 10 m, estimate 
the mass M of the Earth. 

(Solution on page 50 | 


It can be shown that the model of a constant gravitational force is quite a reasonable 
one for heights of a few kilometres above the Earth's surface. Consider two points P; 
and P» lying on the same radial line, with P; on the surface and P, at a distance h 
above it (see Figure 9). We wish to compare the force of gravity at points P; and P2 on 
a particle of mass m. 


At both points the force on m is ‘downward’, along the direction of —e,. The 
magnitude of the force at P, on the Earth’s surface (at radius R) is 


_ mMG 
= as 


It is said in some quarters 
that it took Newton twenty 
years to prove this result, 
causing him to delay 
publishing Principia. But 
this version of events has not 
been established to the 
satisfaction of all science 
historians, 


Figure 9 
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At the point P2, which is a height h above P;, the force has magnitude 


mMG 
oe Cece 
The ratio of these two magnitudes is 
Fy R 1 


FR (R+hP (+h/R? 
Provided that / is small compared to R, the force of gravity at point P» is nearly equal 
to that at P,. 
Exercise 16 


Given that the acceleration due to gravity is 9.81ms~? at the Barth's surface, and that the 
radius of the Earth is 6.4 x 10° m, estimate the acceleration due to gravity at a height of 104m 
above the Earth’s surface. 


[Solution on page 50 | 


Summary of Section 2 


1. Kepler’s laws of planetary motion are as follows. 
(I) Each planet moves in an ellipse, with the Sun at one focus. 
(II) The line joining a planet to the Sun sweeps out equal areas in equal times. 


(III) The square of the orbital period of a planet is proportional to the cube of the 
semi-major axis of its orbit. 


2. (i) The ellipse shown in Figure 10 has semi-major aris a and semi-minor aris b. 
With respect to the axes OXY, its equation is 


Be 
eR 


=1, 


where a > b> 0. 


Figure 10 


(ii) The area of the ellipse is rab. 


(iii) The foci of the ellipse aré located a distance q to either side of the centre O', 
where 


(iv) With respect to axes Oxry with origin at the right-hand focus, and in polar 
coordinates, the equation of the ellipse is 


Fa 1400s, 
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where 
ra 
a 


1 


is the semi-latus rectum, and 


is the eccentricity, where 0 < ¢ <1. 


3. Newton’s universal law of gravitation states that the gravitational force on a 
particle with mass +m, and position vector r; due to another particle with mass m. if 
and position vector ry is 


Gmyum. 
Fy = ——(rj - 11), 
y pene 37%) 


where G is the gravitational constant, 


4, A body has a spherically symmetric mass distribution if (with the origin of 
coordinates at the centre of the body) its mass density function p(r) depends only 
on the radial distance |r| =r = \/x? + y? + 2?. The gravitational force exerted by 
such a body B on an exterior particle with mass m and position vector r is given 
by 


GmM 
7a ers 


F(r) = —gradU(r) = — 


where e, = r/r is the unit vector in the outward radial direction, 
GmM 


U(r) =-=— 


is the gravitational potential energy of the particle, and 
M= vf p(x!) dV 
B 


is the total mass of the body, 


3 (Orbits 


In Subsection 2.2 Newton’s universal law of gravitation was introduced, leading to a 
derivation of the corresponding gravitational potential energy function for a spherically 
symmetric mass distribution. In the current section we shall use this result in 
developing a mathematical model to describe the orbit of a body under the 
gravitational attraction of a much larger body. To a reasonably good approximation, 
this model applies to the orbits of a planet around the Sun, of a moon around a planet, 
of a man-made satellite around a planet or around the Sun, and of comets around the 
Sun. The model will enable us to derive Kepler’s three laws of planetary motion, which 
were stated and explained in Subsection 2.1. 


The modelling assumptions which suffice to yield these laws are as follows. 


Modelling assumptions for gravitational orbits 


(i) The world-view of Newton is valid, including his three laws of mechanics 
and his universal law of gravitation. 


(ii) All forces other than the gravitational attraction between the two bodies 
may be neglected. (In particular, the gravitational effect of other 
celestial objects is ignored.) 


(iii) One body is much more massive that the other. 


(iv) The more massive body has a spherically symmetric distribution of mass. 


Kepler's laws are also the 
first item in the Summary of 
Section 2 on the previous 
page. 
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Of these assumptions, (i) is the hardest to improve upon and (iii) is the easiest, 

while (ii) and (iv) are intermediate. It was only with Hinstein’s general theory of 
relativity in the early twentieth century that small corrections to Newton’s world-view 
(Assumption (i)) were introduced. This was done in order to explain small departures 
from the Newtonian prediction that had been observed in the orbit of the planet 
Mercury. We could (but will not!) improve upon Assumption (iii) using the 
mathematics developed in the course. This improvement is required, for example, to 
describe the orbits of two stars of comparable mass under their mutual gravitational 
attraction. 


3.1 Central forces 


By Assumption (iii) above, we suppose that the more massive body (of mass M) does 
not move at all under the gravitational action of the lighter body (of mass m). Its 
centre of mass coincides with its geometric centre, by Assumption (iv). For convenience, 
we choose the origin of coordinates to be at this centre, as shown in Figure 1. 


Figure 1 


It was shown in Unit 17 that the centre of mass of a system of particles acted upon by 
external forces moves as if all the mass of the system were concentrated there. We may 
therefore treat the motion of the centre of mass of the lighter body as if it were a 
particle. Furthermore, since the heavier body is assumed to be spherically symmetric 
(in gravitational terms), we can apply the result of Subsection 2.2 which equates the 
gravitational effect of such a body with that of a particle of mass M at its centre. The 
resulting two-particle view is shown in Figure 2, where r is the position vector of the 
lighter particle. + 


There is a simplifying feature of this model which can be shown at the outset: the 
motion of m occurs in a plane. The proof of this statement rests on the fact that the 
gravitational force exerted on m is a central force. In general, this term describes a 
force which is directed from the particle m towards or away from the source of the 
force. If the source is at the origin, this means that a central force has the direction of 
—r or of r, The definition is as follows. 


A central force F, acting from the origin on a particle with position vector 
r, has the form 


F(r) = g(r), () 


where g(r) = 9(«,y,2) is a scalar function of position. 


Exercise 1 


Which of the following are central forces (where K is a constant)? 


(@) F(@)=Kr (@) Fe) = “45 
(b) F()=4r @ Fe)=KG+i+K 
© Fe=“r (9) Fe) = “(i+ w+) 


Unit 17 Subsection 2,3 


A similar argument to this 
was employed in the solution 
to Exercise 14 in Section 2, 


Equation (19) of Section 2. 


As earlier in the unit, we use 
the symbol m as a convenient 
label for the particle of mass 
m. 


fA 
Figure 2 
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Exercise 2 
What is the scalar function g(r) for the gravitational system of Figure 2? 
[Solutions on page 50 | 


The angular momentum L of a particle with mass rm and position vector r relative to 
the origin is, as defined in Unit 29, 


L=rx mi. (2) 


This definition is not confined to two-dimensional motion, though the main application 
in Unit 28 was to motion in a circle, However, it is a consequence of the results of 
Exercise 3 below that if the particle is acted upon by a central force then its motion 
does lie in a plane. 


Exercise 3 
Show that: 


(i) when the angujar momentum L is non-zero, the position vector r is always perpendicular 
to L; 


(ii) when the force on the particle of mass m is a central force then L = 0, so that L is a 
constant vector; 


(iii) when L = 0, the velocity # and position vector r are parallel or anti-parallel if neither is 


20r0. 


[Solution on page 50] 


From parts (i) and (ii) of Exercise 3, we can say that if L 4 0 then the position vector r 
of the particle of mass m is always perpendicular to the fixed direction defined by the 
constant angular momentum vector L. In the general case then, when L is non-zero, 
the particle acted upon by a central force moves in the plane perpendicular to the fixed 
direction of L, as shown in Figure 3. Also (from part (iii)), in the event that L is zero, 
the velocity lies along a straight line from the origin, and the particle will continue to 
move along this straight line. 


Figure 3 


The constancy of the angulag momentum L for a particle in gravitational orbit arises 
from the fact that gravity is a central force, and not from any other detail of its form. 
This result applies also, therefore, to other central forces. 


Motion under a central force 


A particle with mass m and position vector r, which moves solely under the 
action of a central force, has constant angular momentum 


L=rx mer. (2) 


Its motion is confined to the plane perpendicular to the direction of L when L 
is non-zero and to a straight line if L = 0. 


Unit 29 Subsection 1.2 


In Unit 29 we used the 
symbol l for the angular 
momentum of a single 
particle and L for the angular 
momentum of a system of 
particles. Throughout this 
unit we refer only to the 
angular momentum of a 
single particle, and for visual 
convenience this is 
represented by the symbol L. 
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As the motion is planar, we can now choose to orient the coordinate axes so that the 

motion lies in the (x, y)-plane, as indicated in Figure 4. This will permit us to apply 

the apparatus of plane polar coordinates which was developed in Subsection 1.3. Tt has yet to be established 
that the planar motion of a 
particle in gravitational orbit 
is as simple as that depicted 
here. For example, some 
other central forces lead to 
spiralling motions. 


uv 


at 


Figure 4 


Since the angular momentum L is a constant vector, both of its direction and 
magnitude are constant. The constancy of its direction sufficed above to show that the 
orbit always lies in a plane. This is actually part of the content of Kepler’s first law, 
which asserts that the orbits are ellipses and therefore plane figures. 


We have not yet used the additional information that the magnitude of L is constant, 
and this is actually sufficient to prove Kepler’s second law. As you may recall from 
Subsection 2.1, this law relates to the rate at which area A is swept out by the line 
joining a planet to the Sun (which corresponds here to the line of the position vector r). 
We showed that A was related to the radial distance r and the angular velocity @ by 
the equation 


A= 4r76. Equation (9) of Section 2 
Kepler's second law states that A is constant, which from above is equivalent to the 
condition 
da, 
= =0. 
F(r*6) @) 
In the exercise below you are asked to show that Kepler’s second law is a consequence 
of the constancy of angular momentum. 


Exercise 4 


(i) By expressing r and # in terms of the plane polar unit vectors e, and eg (Equations (8) 
and (12) of Section 1), show that the angular momentum of a particle moving in the 


(x, y)-plane is L = L.k, where 


L. = mr". 


(ii) Deduce that if Z. = 0 then Equation (3) is satisfied. 
[Solution on page 50] 


In summary, thus far, we can say the following. 


Kepler's second law, and the fact that each planet moves in a plane, are 
consequences only of the fact that the gravitational force between the planet 
and the Sun is a central force; taken together, they amount to a statement of 
the constancy of angular momentum. 
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3.2 Isotropic central forces and gravity 


With the force of gravity in mind as the ultimate target, we now specialize from central 
forces to the more restricted case of isotropic central forces. These are central forces for 
which the magnitude depends only on the radial distance from the origin. 


An isotropic central force has the form 
F(r) = f(r)er. (4) 
where r = \/x? + y? + 2* is the distance from the origin, and e, = r/r. 


Exercise 5 
(i) Which of the forces listed in Exercise 1 (page 28) are isotropic central forces? 


(ii) Show that the gravitational force is isotropic, and find the function f(r) for this case. 


Solution on page 50 | 


Choosing our axes as in Figure 4, so that the orbit lies in the (7, y)-plane, we have 


0 and r= \/x? + y?. Newton's second law, for the motion of a particle of mass m 
acted upon by an isotropic central force, is 


#= f(re,. 


Since the right-hand side of|this equation is directed along the unit vector e,, it makes 
sense to express the acceleration # in terms of plane polar coordinates, using 

Equation (13) of Section 1. The particle's equation of motion in the plane can then be 
written as 


m 


‘ d fi 
m(i = ré)e, + alr Deo =f(r)er. 


from which we have 


m(# — 0?) = f(r) (5) 
md on _ 
and. road 0) =0. 


The second equation is equivalent to Equation (3), whose validity was established for 
any central force in Exercise 4, The consequent constancy of the angular momentum 
magnitude, 


L. =mr*6, (6) 


can be exploited to eliminate the variable @ in Equation (5). Thus on putting 
4 = L./mr? into Equation (5), we obtain 


mi — mer ( 


or mi = f(r) + 


(7) 


Equations (6) and (7) contain, between them, the possibility of a complete solution to 
the problem of finding the orbit of a particle of mass m under the action of any isotropic 
central force, given by Equation (4). The procedure to do this would be as follows. 


mr 


(i) The initial conditions are the values of [r,4] and {i,6] at some initial time, say 
t=0. 


(ii) Since L. is a constant, these initial conditions determine the value of L. for all 
subsequent times, from Equation (6). 


(iii) Equation (7) may be solved, either analytically if possible or numerically, to obtain 
an expression or values for r(t), where t > 0. 


(iv) This solution may be substituted into Equation (6) to obtain an expresion for 0(t), 
from which 6(t) can be found by direct integration. 


A central force was defined 
by F(r) = 9(r)r, which is 
equivalent to 

F(r) = g(r)rey. 
This force is also isotropic if 
g(r) = g(r). 
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The approach we shall eventually adopt, however, is somewhat different. Instead of 
seeking r(t) and 0(t) individually, which are parametric equations for the orbit, we look 
for a direct relationship r = r(@). This may describe an orbit which is bound or 
unbound. For a bound orbit (see curve (a) in Figure 5), the distance of the particle 
from the origin O is bounded for all time. For an unbound orbit (see curve (b) in 
Figure 5), the particle approaches from infinitely far away, is deflected by the force 
centred at O and departs again to infinite distances. Corresponding to the bound orbit 
shown in Figure 6(a), the function r(@) has the graph of Figure 6(b). 


Figure 6 


Thus far, the results obtained relate to motion under any isotropic central force, and 
not just to gravity. At this point we shall finally specialize to the case of gravitation, 
though the results may also be applied to any other force which is described by a 
similar inverse square law. 


Specializing to gravity means using the fact that for this particular isotropic central 
force we have (from Exercise 5(ii)) 


GmM 
fein See, 


so that Equation (7) for the particle motion becomes 
2 GmM 


mit = AE, — SOM (8) 
This equation looks somewhat formidable, but it can be tamed, In Subsection 3.3 you 
will see how a relationship r = r(@) can be obtained from Equation (8) to describe 
gravitational orbits. 


For the remainder of this subsection we shall investigate what can be deduced about 
the motion without solving Equation (8) directly. This indirect approach depends on 
the constancy of both the angular momentum L. and the total mechanical energy, EL. 


It follows from the result of Exercise 13 of Section 1 that, for a particle moving in two 
dimensions and acted upon by the (conservative) gravitational force, the constant total 
mechanical energy Z may be written as 


EB = 4n(i* +176) + U(r), (9) 


where U(r) = —GmM/r is the gravitational potential energy function. From 
Equation (6), the constant angular momentum is given by L, = mr?0. After using this 
equation to replace @ by L./(mr?), Equation (9) becomes 


3 _ GmM 


B=3 


(10) 


Hence, by using the constancy of the angular momentum, we have expressed the total 
mechanical energy in terms of 7 and r only. 


fa) {b) 
Figure 5 


The electrostatic force is one 
example of an isotropic 
central force. You considered 
a two-dimensional example of 
such a force in Exercise 15 of 
Section 1. 


Equation (19) of Section 2 


We write U(r) rather than 
U(r) here to emphasize that 
the gravitational potential 
energy U depends only upon 
the radial distance r. 


Tn fact, any isotropic central 
force is conservative, and has 
a constant total mechanical 
energy E given by 

Equation (9). You will see 
this shown in Section 4. 
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As a result, Equation (10) has the appearance of an energy equation for a particle 
moving in one dimension. The first term, }mi?, resembles the one-dimensional kinetic 
energy. Since the bracketed term in Equation (10) looks like a one-dimensional 
potential energy function, we shall call it the effective potential energy and denote it by 
U**(r), This gives 


E = mi? +" (r), 
where the effective potential energy is 
2 GmM 


2mr? r 
In summary, then, the radial motion of the particle is as if it were moving in one 

dimension, with r > 0, under the influence of a one-dimensional conservative force 
having potential energy function U*"(r). From this insight you might expect that 
Equation (8) for the radial motion could be written as 


u(r) = 


(11) 


d 
f= —-—U"(p), 
mi = - u(r) 


by analogy with the definition of a one-dimensional potential energy function. This is 
indeed the case, as can be verified by differentiating Equation (11). 


Obviously, the particle whose orbit we seek is not really moving in a single dimension, 
since as 7 varies with time the angle @ is changing too. However, we have managed in 
Equation (10) to express the total mechanical energy of the orbiting particle in the 
same format as that for a one-dimensional motion, This means that we can deduce 
information about the radial coordinate of the orbit, using the effective potential 
energy function, with the approach that was adopted in Subsection 1.1 for an actual 
one-dimensional motion, 


You saw in the analysis of Subsection 1.1 that whether a motion is bound or not 
depends both on the shape of the graph of the potential energy function and upon the 
value of the constant E. Here the effective potential energy is given by Equation (11), 
and has the graph shown in Figure 7. The main features of the function U' illustrated 
here are as follows. 


(i) When r is small, the term L?/(2mr?) dominates the term —GmM/r, so U""(r) 
becomes large and positive as r approaches zero. 


(ii) When r is large, both of the terms in U*"(r) tend towards zero, so that the limit of 


u*"(r) for large values of r is zero. 


(iii) For a range of intermediate values of r, the term —GmM/r dominates the term 
L2/(2mr?), causing U*"(r) to be negative with a minimum value at some point 
r=r, 


Exercise 6 


Find the value r =r" for which U*™ has a minimum, and calculate the corresponding 
minimum value UST, = U(r), 


[Solution on page 50] 


Figure 8 shows additions to the graph of U*" relating to two values of the total 
mechanical energy £, which correspond to a bound and an unbound motion. A particle 
approaching the origin O with some positive energy, E""°""4, reaches a closest point 
to O at r =a, following which it is always moving away towards larger values of r. On 
the other hand, a particle with some negative energy, B?°'"4, undergoes a periodic 
oscillation of the radial coordinate between a minimum value r = b and a maximum 
value r = c. (In Subsection 1.1 these would have been bound and unbound motions 
along a line. Here they correspond to bound and unbound orbits.) 


Exercise 7 


Describe the orbit corresponding to the total mechanical energy B = US 


min* 


[Solution on page 50 | 


‘As a check, if Equation (10) 
is differentiated with respect 
to time (with B = 0), then 
the outcome is Equation (8). 
Put another way, 

Equation (10) is a first 
integral of Equation (8). 


Equation (1) of Section 1 


ye 


ve, 
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The following summary details what information has been obtained so far about orbital 
motion under the influence of gravity. 


Gravitational orbits 


Consider the planar gravitational orbit of a particle, with mass m and polar 
coordinates [r, 4], due to a spherically symmetric body of much larger mass 
M, situated at the origin. 


(i) The angular momentum L. of the particle is constant, where 
L.=mr*6. (6) 
(This is a property of any central force.) 
(ii) The total mechanical energy E of the particle is constant. It is given by 


(a) E = jm(r? +176?) + U(r), (9) 
where U(r) = —GmM_/r is the (actual) potential energy of the 
particle; 

(b) E = }mi?* +U%(r), 
where 

uet(r) = me Ss iu (11) 


is the effective potential energy. 


(iii) The orbit is unbound if £ > 0, and bound if Us, < E <0, where Us, 
is the minimum value of UeH(y), 


(iv) The equation of motion for r(t) is 


L? | GmM 
S- 


m= 


é (8) 


— 


The properties of gravitational orbits arising from the constancy of angular momentum 
and of total mechanical energy can be applied to obtain useful information in specific 
cases, as illustrated by the exercises below. 


Exercise 8 


A particle of mass m is launched with speed v at right angles to the line which joins it to a 
spherically symmetric body of much larger mass M. At the instant of launch, the distance 
between m and the centre of M is R. 


(i) Find the angular momentum of the particle for all subsequent times. 
(ii) Find the total mechanical energy of the particle for all subsequent times, 
(iii) Show that the motion will be bound provided that v? <2MG/R. 


Exercise 9 - 

A particle of mass m is launched as described in Exercise 8. Its speed of launch y and distance 
of launch R are related by the condition v? = MG/(2R), 

(i) Find the total mechanical energy, and show that the resulting orbit is bound. 


(ii) Find rmin and rmax, the respective smallest and largest distances of the particle from the 
centre of the attracting mass M. (Hint: From the solution to Exercise 8, the total 
mechanical energy is E = mv? ~ GmM/R, but this is also expressible as 
E = mi? + U(r) for all values of r and # on the particle's trajectory. At the turning 
points, when r = rin or r = Tmax, we have 7 = 0. Hence the equation for the turning 
points is F = UM (p).) 


[Solutions on page 50 | 


Here rin and rmax 
correspond to the values b 
and ¢ of Figure 8. 
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3.3 The equation of gravitational orbits 


We now revert to consideration of the radial equation of gravitational motion, 


. 22  GmM 
i = (8) 


As mentioned earlier, we seek to obtain from this equation a relation of the form 

r =1(6) to describe gravitational orbits. To do so, we apply the chain rule of 
differentiation to convert Equation (8), which specifies a function r(t), to an equation 
for r(@). For the first derivative, #, we have 

dr _ ddr _ 5 
dt di dd do” 

From Equation (6), 6 can be replaced by L./(mr?), where L. is a constant. The 
previous equation can then be expressed as 


L. dr 


"mr? dd 


or, in a form which will shortly be convenient, as 


1 
mae (7): (2) 


The equation of motion contains 7, which may be written in terms of a derivative with 
respect to @ by a further application of the chain rule. Thus, we obtain 


d,_d0d._ id. 
~ dt’ edb’ do 
which, after another use of the relation 6 = L./(mr*) and substitution for * from 
Equation (12), becomes 


se () 


“mtr 
Equation (8) can now be written as 


Defining 


u(0) = ea (13) 


gives, after some rearrangement, 


Pu 


ae t u(@) = —— 


(14) 


Equations (13) and (14) contain the information necessary to deduce the orbit r = r(@). 
You have met equations like Equation (14) many times. It is a second-order, 
constant-coefficient, inhomogeneous differential equation, whose general solution has 
two arbitrary constants. One way of writing the solution is 


m?MG 
be” 
‘The two arbitrary constants in this solution are the amplitude B and the phase @. (We 
assume that B is non-negative. This represents no loss of generality, since the sign of 
the term B cos(@ — @) may be changed by adding 7 to the arbitrary phase @.) The 


orbital equation is obtained from Equation (15) on using Equation (13) to replace u by 
1/r, After multiplying through by 


u(0) = Bcos(@ — 09) + (15) 


LB 


(16) 


and putting e = Bl, the orbital equation can be written as 
U 


Fm = 1 te c0s(0 ~ H). (17) 


For example, in Unit 7 
Subsection 2.4 a similar 
equation was derived when 
modelling the motion of a 
particle suspended from a 
perfect spring in a uniform 
gravitational field, 


Here B is a non-negative but 
otherwise arbitrary constant, 
so the same is true of . 
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Equations (16) and (17), together with the knowledge that L- (or, equivalently, 1) and 
e are constants, provide the required mathematical description of gravitational orbits 
for the model outlined at the beginning of this section. The remainder of the section 
will be spent in examining the implications of these results. 


Note first that the introduction of the symbols / and ¢ has permitted us to write the 
orbital equation in a form which is identical to Equation (8) of Section 2. Provided that 
e <1, therefore, Equation (17) represents an ellipse with eccentricity e and semi-latus The terms ‘eccentricity’ and 


rectum 1. ‘semi-latus rectum’ are used 
to describe e and | whether 
Exercise 10 or not the orbit is an ellipse. 


Equation (17) describes a general gravitational orbit. 


(i) By considering the possible values of the right-hand side of this equation, show that the 
distance of closest. approach of the particle of mass m to the origin occurs at the angle 
4 = %, and has the value 


U 
T¥e’ 
(ii) Show that, when e < 1, the maximum distance of the particle of mass m from the origin 
occurs at the angle @ = +7, and has the value 
a 
l-e 
(iii) Show that when e > 1 there is no maximum value of r. 


Tin = 


Tinax.= 


(iv) From parts (i)-(iii), give conditions on the eccentricity e which characterize bound and 
unbound orbits, 


[Solution on page 51] 


From the results of Exercise 10 we can sketch the possible types of orbit that 
Equation (17) represents. This is done in Figure 9 below. When e < 1, the orbits are 
closed (ellipses) with closest and furthest distances ryyin and ryax respectively. When 
€ 2 1, the orbits are unbound (hyperbolas or parabolas), and only the distance of 
closest approach ryjn is defined. In either case, the orbit is symmetrical about the line 
0=6. 


(ajo<d (b)e >a 


Figure 9 


It is convenient to consider the gravitational orbits with respect to a set of axes which 
are rotated anticlockwise by an angle 0 from those of Figure 9, This is equivalent to 
choosing % = 0 in the orbital equation, or to choosing the x-axis to coincide with an 
orbital symmetry, and represents no loss of generality. In place of Equation (17) we 
then have 


t 
7 =i +ecosd, (18) 


where the semi-latus rectum / is given as before by Equation (16), The corresponding 
pictures for the orbits are shown in Figure 10. 
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Figure 10 


Equation of gravitational orbits 


The Newtonian model for the gravitational orbit of a body of mass m around 
a spherically symmetric body of much larger mass M located at the origin 
leads to the orbital equation 


t 
= = 1+ ecos6, (18) 
where the eccentricity ¢ is non-negative and the semi-latus rectum / is given 
by 
2 
l=—. 16 
m?MG (16) 


Here L. is the constant angular momentum magnitude of the orbiting body. 
If 0 < e <1 then Equation (18) describes an ellipse. 


Values for the constants / and ¢ in Equation (18) can be calculated in a Specific case 
from experimental observations, and this can be done in various ways. For example, 
suppose that the distance of closest approach rj, and the corresponding angular 
velocity @ have been measured, and that the value of the central mass M is known. 
Then a value for L./m can be obtained from Equation (6), This leads in turn to a 
value for | from Equation (16) and a value for e from the result of Exercise 10(i). 


The energy analysis carried out in Subsection 3.2 predicted that gravitational motion 
will be bound if the total mechanical energy satisfies the inequality US", < E <0, 


min S 
and unbound if E > 0, It is to be expected that these predictions are related to those 


obtained in terms of the eccentricity e in Exercise 10, namely, that an orbit is bound if 
e <1 and unbound if e > 1. We consider this relationship in the example below. 
Example 1 


Show that the total mechanical energy can be written as 


Hence verify that the condition E > 0 for unbound orbits is equivalent to e > 1, and 
that the condition E < 0 for bound orbits is equivalent to e < 1 (given that e > 0). 


Solution 


From Equation (10), 2 can be written as 


This analysis was based on 
Equations (10) and (11), 
together with Figure 8. 
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Using Equation (12) for 7, this becomes 
2[d/1\]?, 12  GmM 
= 1m— |—[-— fain 2) 
ti E (*)] Tom? 


From Equation (16) for /, this is 


e982 ([0 +0) -O} 


But //r is given in terms of @ by Equation (18), so that 


B= PMG (62 sin® 9 + (1-4 ecos0)? ~ 2(1 + ¢0086)}. 
On simplifying, this gives 
mMG_ 4 
=3@-0 


So the condition E > 0 for unbound orbits corresponds to e > 1, and the condition 
E <0 for bound orbits corresponds to e < 1. (Although this outcome is expected, it 
represents a useful cross-check on the results achieved.) © 


Exercise 11 

Describe the orbit given by Equation (18) when the eccentricity ¢ is zero, What are rmax and 
Tmin in this case, and what is the total mechanical energy? 

Exercise 12 


A satellite of mass m follows a circular orbit with constant angular velocity and radius d 
around a spherically symmetric planet having a much larger mass M. 


(i) Show, from first principles (using Equation (13) of Section 1), that the satellite's 
acceleration lies along the radius of its orbit and has magnitude 


le] = a6. 
(ii) Show too that its orbital angular momentum has magnitude 
L, = md*6. 


(iii) By equating mi to the gravitational attraction on the satellite (Newton's second law) 
show that the value of the orbital radius agrees with that found in Exercise 11, 


(iv) By adding the kinetic energy of the satellite to its gravitational potential energy, verify 
that the value of the total mechanical energy agrees with that found in Exercise 11. 


[Solutions on page 51 | 


3.4 Planets and comets 


In this brief subsection we shall verify the statement made at the beginning of this 
section to the effect that Kepler's laws of planetary motion are consequences of our 
model based on Newton’s laws. Following this there is a short discussion about the 
orbits of comets. 


Kepler's first law says that each planetary orbit is an ellipse, with the Sun at one focus. 
You showed in Exercise 3 that a planet’s orbit lies in a plane because gravity is a 
central force. According to the model, a gravitational orbit within this plane has an 
equation of the form 

l 

5 =i tecosd, (18) 
where the semi-latus rectum / and the eccentricity e are positive constants for a 
particular planet. 


Planetary orbits are clearly bound, and you established in Exercise 10 that e <1 is a 
necessary and sufficient condition for bound motion. But for e < 1, Equation (18) is the 
equation developed in Subsection 2.1 to describe an ellipse with one focus at the origin. 
In our Newtonian model of planetary motion, the Sun is located at this focus. Thus 
Kepler's first law is a consequence of the model. 


Equation (7) of Section 2 
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Kepler's second law says that the line joining a planet to the Sun sweeps out equal areas 
in equal times. We showed in Subsection 3.1 that this is equivalent to the equation 


£8) =0, (3) 


which expresses the constancy of the planet’s angular momentum. This is again a 
consequence of gravity being a central force. 


Kepler's third law is a relation between the period of a (bound) orbit and its lateral 
extent. It says that if T is the period and a is the semi-major axis of the elliptical orbit, 
then 


T? = ka’, (19) 


where k is a constant not depending on T or on a. In Exercise 13 below you are asked 
to show from basic principles that Equation (19) holds for the special case of a circular 
orbit, where the eccentricity e is zero. In this case the semi-major and semi-minor axes 
both equal the radius of the circle. In Exercise 14 you are then asked to establish 
Kepler's third law for any elliptical orbit. 


Exercise 13 
(i) Show that the magnitude of the gravitational force between the Moon (of mass m) and 
the Earth (of mass M) is 
GmM 
a 
where d is the separation between their centres. 


(ii) Assume that the Moon is moving uniformly in a circle around the Earth, so that (from 
Exercise 12(i)) the magnitude of its acceleration is || = d6*, where 4 is its constant 
angular velocity. Show that 


le] = An?d 


where T is the period of the Moon's orbit. 


(iil) Show that if air resistance is ignored then the acceleration of an apple dropped close to 
the surface of the Barth is 


where R is the Earth's radius. 


(iv) From parts (i)-(iii) show that, according to the Newtonian model of gravitation, 

Ant 

gh 

(vy) It was known by Newton that d ~ 60R, g ~ 9.8ms~? and T' ~ 28 days. but it is now 
believed that he used for R, the Earth's radius, a value of 3500 ‘Italian miles’, where one 
Italian mile is 5000 feet (rather than 5280 feet), This would make R ~ 5.3 = 10° m. Using 
these values (converted to appropriate SI units where necessary), do you find that the 
left- and right-hand sides of the predicted relation in part (iv) agree ‘pretty nearly’? (A 
more accurate value for R is 6.4 x 10° m.) 


T* =kd*, where k 


Exercise 14 


Consider an elliptical orbit with semi-major axis a and semi-minor axis b. Then it is known 

from Exercise 2 of Section 2 that the area of the ellipse is tab. Also, the semi-latus rectum is 
a nMG 

(from Equation (5) of Section 2 and Equation (16) of this section), and the rate at which area 

is swept out is 


A=tr6 


L. 

2m 

(from Equation (9) of Section 2 and Equation (6) of this section). Here L. is the (constant) 
orbital angular momentum. 


(i) From these facts, show that the orbital period, 7’, and the semi-major axis, a, are related 
by Kepler's third law, in the form 


An® 
(ii) Show that, for a circular orbit around the Earth, this agrees with the result of 
Exercise 13(iv). 


T 


ka*, where k 


[Solutions on page 51 | 


This exercise refers to 
Newton's quotation given at 
the beginning of 

Subsection 2.2. 


40 MST 204 30,3 


Kepler's laws were formulated to describe the bound orbits of planets, so there is no 
reason to suppose that they should have significance for unbound orbits. The first and 
third laws can by their very nature apply only to the bound case. However, Kepler’s 
second law applies even for unbound orbits since, as you have seen, it depends only 
upon the fact that gravity is a central force. 


The Newtonian model which led to Equation (18) for gravitational orbits makes no 
assumption as to whether or not the motion is bound (with eccentricity e < 1) or 
unbound (with e > 1), As shown in Example 1, these conditions for the eccentricity are 
related to similar conditions on the total mechanical energy E of an orbit. If E < 0 for 
a body moving under the gravitational influence of the Sun, then the body has 
insufficient energy to ‘escape’ from the Sun's gravitational field and is confined to 
bound and periodic motion. Bodies with energy Z > 0 may approach the Sun once, but 
are not ‘captured’ by it. These unbound orbits are called parabolic when e = 1 (E = 0) 
and hyperbolic when e > 1 (E > 0). 


At any particular time there are hundreds of bodies in the solar system moving in nearly 
or actually unbound orbits under the Sun’s gravity. These are the comets. Precise 
details of cometary orbits are not easily obtained, but it appears that most comets are 
in weakly bound orbits, with values of eccentricity ¢ just less than unity. These comets 
therefore have highly eccentric, or flattened, orbits, as illustrated in Figure 11. 


u 
b 
4 - 
—<$———ae 7) 7° 


As a first approximation, the orbits of these bound comets obey Kepler's laws, but as 
they follow their elongated trajectories through the solar system they occasionally pass 
close to a planet. In these circumstances the gravitational effect of this third body can 
alter their orbits, even to the extent of raising their energies to positive values. Owing 
to its relatively large mass, Jupiter is a particular offender in this regard. The most 
famous bound comet is Halley’s comet, which has a period of about 76 years and last 
passed close to the Sun in early 1986. Its orbit has eccentricity e ~ 0.97 and a 
semi-major axis of about 18 times the distance from the Barth to the Sun. Comparable 
data for the planets, the Moon and Pluto are given in the table below. 


Figure 11 


Name Mass Period of Semi-major Eccentricity Mean 
revolution axis of orbit of orbit diameter 
around Sun 
Sun 332:948.0 109.1 
Moon 0.012 0,055 0.271 
Mercury 0.055 0.241 0.387 0,206 0.391 
Venus 0.815 0.615 0.723 0,007 0.951 
Earth 5.974 x 1074 kp | 365.256 days | 149.6 x 10°km 0.017 12735 km 
Mars 0.107 1.88 1.524 0.093 0,531 
Jupiter 317.9 11.86 5.203 0.048 10.98 
Saturn 95.1 29.46 9.539 0.056 9.09 
Uranus 14.6 84.00 19.182 0.047 4.00 
Neptune 17.2 164.79 30.058 0.009 3.78 
Pluto 0.17 247.70 39.44 0.254 0.22 


At any time there are comets in the solar system that are pursuing unbound 
trajectories. These bodies come in from distant places, pass partly around the Sun to a 
distance of closest approach and then recede, never to return to our solar system. 


You investigated the 
connection between the 
eccentricity and shape of an 
ellipse in Exercise 4 of 
Section 2. 


For each body except the 
Earth, the values given in the 
table for mass, period, 
semi-major axis and mean 
diameter are multiples of the 
corresponding values for the 
Earth. 
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Summary of Section 3 


iS 


A central force F, acting from the origin on a particle with position vector r, has 
the form 


F(r) = o(t)r, 


where g(r) = g(x,y, 2) is a scalar function of position. If the particle has mass m 
and is not acted upon by any other force, then it has constant angular momentum 
L=r x mg. If the unit vector k is chosen in the direction of L, then 


(i) the motion of the particle is confined to the (2, y)-plane; 
(ii) L=L.k, where L. = mr?6 is constant. 


The gravitational force due to a spherically symmetric body whose centre is fixed 
at the origin is a central force. In this case, the constancy of L, is equivalent to 
Kepler’s second law. 


An isotropic central force has the form 
F(r) = f(r)er, 


where r= \/2? + y* + 2? is the distance from the origin, and e, = r/r. For a 
particle of mass m acted upon by such a force, the radial equation of motion is 


s i 
mi = f(r) + = 
which may be integrated (analytically or numerically) to give r(t). The 
gravitational force is isotropic, with f(r) = -GmM/r?, where M is the central 


mass and G is the gravitational constant. (It is assumed that M is much larger 
than m, in order that the central mass should remain fixed at the origin.) 


For the particle moving under gravity as described above, the total mechanical 
energy F is constant. It may be expressed as 


(i) E = }m(i +176") + U(r), 
where U(r) = —GmM/r is the (actual) potential energy of the particle; 
(ii) B= }mi? +U"(r), 
where 
12 GmM 
petty.) ee 
OSG bear tar 


is the effective potential energy. 


Orbits are unbound if £ > 0, and bound if U«", < 2 < 0, where US, is the 
minimum value of U(r). 
The equation of gravitational orbit for the particle described above is 
tf 
-~=1+ecos0, 
r 
where the semi-latus rectum !/is given by 
es 
m?MG" 
If 0 <e <1 then the orbit is an ellipse with eccentricity e, showing that Kepler's 
first law is satisfied. If e > 1 then the motion is unbound. 


In terms of the eccentricity and the semi-latus rectum, the total mechanical energy 
is 


_ GmM 
~~ Bf 


(e? -1). 


Kepler’s third law is also a consequence of this Newtonian model, since the period 
T of an elliptical orbit and its semi-major axis a are related by the equation 

4x? 

MG 


3 


a where k= 
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4 Central forces and conservation laws 
(Television Section) 


4.1 Isotropic central forces 


‘The television programme considers some general statements that can be made about 
the motion of a particle in three dimensions which is acted upon by a central force and, 
more especially, by an isotropic central force. When the force is central, then the 
angular momentum of the particle is constant and its motion lies in a plane. When in 
addition the force is isotropic, one can say further that there is a potential energy 
function for the particle, so that it has a constant total mechanical energy. 


The first of these results was established in Subsection 3.1. The second was applied for 
the particular case of gravity in Subsection 3.2, but holds more generally. This more 
general version is a consequence of Result (i) below, 


(i) Ifa scalar function @(r) depends upon position only through the radial 
distance |r| = r = (a? +y* + 2*)!/2, so that #(r) = (r), then the 
gradient of @ can be written as 


grad = oe A 
dr 


where e, = r/r is the unit vector in the outward radial direction. 
(ii) Hence if F is an isotropic central force, that is 
F(r) = f(r)e, 


for some function f, then F is conservative, with a potential energy 
function U(r) such that 


The motion of a particle in this force field is planar, and has constant 
total mechanical energy 


B= 4m(i? + 76?) + U(r). 


Result (ii) follows from Result (i) on referring to the definition of a conservative force 
(Equation (3) of Section 1) and to the formula for E derived in Exercise 13 of 
Section 1, (The motion is planar because the force is central, as in Subsection 3.1.) 


The proof of Result (i) is achieved by applying the definition of the gradient and the 
chain rule. Thus we have 


oO, Ly do 


grad @ = 


Since #(r) = #(r), this can be written as 
do dr, dbdr, , dbOr, db (dr, Or, Or 
dr Oz dr Oy)? dr de. dr (G+ dy) * de k) 
But from r = (x? + y? + =*)!/2, we obtain 
ar 
Similarly, we find that 


grad @ = 


ne P 2 2 2 2)-1/2 
x (a? + y? + 2? an? + y? $22)? x on 


dor dé 


dee 
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This completes the proof of Result (i). 


Now watch the television programme. 


4.2 Programme synopsis 


Tom Smith introduces central forces, As an example, he demonstrates an almost 
frictionless puck, connected by a spring to a fixed origin but otherwise free to move on 
a horizontal surface. A central force acting on a particle is a force which is directed 
either parallel or antiparallel to its position vector r. Such a force can be written 
mathematically as 


F(r) = g(r)r, 


where g(r) = g(x.y,=) is a function of position. It is shown that the angular 
momentum of the particle, L = r x mr, is conserved (constant in time) if the particle is 
acted upon only by a central force. In this case the particle moves in a plane, which 
may be taken as the (x, y)-plane. With this choice of axes, the const.nt vector L is 
directed along the z-axis. The fact that a force is central thus allows something to be 
said about the corresponding particle motion, even when it is not possible to solve its 
equations of motion analytically. 


Mike Crampin takes up the story by pointing out that, since L is constant (L = 0), 
both its direction and magnitude are constant. The motion is in the (rc, y)-plane and, 
using plane polar coordinates, the angular momentum is 


L=L.k, where L.=mr*d. 


Since the magnitude L. is constant, this relation can be used to eliminate either r or 6 
from the equations of motion. The constancy of r20 is equivalent to Kepler's second 
law, which says (when the force involved is that of gravity) that equal areas are swept 
out in equal times. So this law of Kepler, together with that part of his first law which 
implies that planetary orbits are planar, follows from the fact that gravity is a central 
force, 


At the next stage, isotropic central forces are introduced, These are central forces 
whose magnitude depends only upon the distance |r| from the origin. Such forces have 
the form 


F(r) = f(r)er. 


where r = |r| and e, is the unit outward radial vector. This is a special case of the 
general central force, g(r)r. 


At this point the second result. proved in Subsection 4.1 is applied, This says that if a 
force is isotropic as well as central then the total mechanical energy of the particle is 
conserved as well as its angular momentum, for then there exists a potential energy 
function U(r) such that 


dU 
=—gradU = ——e,. 
F grad U apt 
Using plane polar coordinates to describe the motion, the total mechanical energy is 
E = 4in(i +16?) + U(r). 


But L. = m4 is constant too, so E can be written as 
L? 
+ un) ; 


2mr? 
This equation expresses the energy solely in terms of r and ?, resembling the energy 
equation for a one-dimensional motion, so the radial motion can be considered 
separately. To make use of the one-dimensional analogy, it is convenient to define an 
effective potential energy function, U", by 


= Ley? 
B= jit + ( 


Then the constant total mechanical energy is 
EB = km? +U"(r), 


Equation (1) of Section 3 


See Subsection 3.1 for details 
of this analysis. It is assumed 
here that L 40. You 
investigated the case L = 0 
in Exercise 3(iii) of Section 3. 


See Exercise 4 of Section 3. 


Equation (4) of Section 3 


In fact, the potential energy 
function is 


u(r) =- f (s)ds, 
Jr 
where ro is any non-negative 
constant. 


The constancy of E (E = 0) 
was shown in Subsection 1.2, 
and its expression in plane 
polar coordinates was derived 
in Exercise 13 of Section 1. 
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Looking at the graph of the function U*", and using the constancy of B, allows us to 
see whether there are any turning points for the radial motion. At these points + 
vanishes, and the radial motion reverses direction. The details depend upon the precise 
form of the potential energy function U and on the value of E. If there are two turning 
points then these are the closest and furthest distances from the origin reached by the 
particle at the given value of energy 2. If there is one turning point only then the 
motion is unbound, and the turning point is the distance of closest approach, 


If the potential energy looks like that in Figure I(a), then U*" is given by the graph in The potential energy given in 
Figure 1(b) and there are always two turning points, for any energy E. But for a Figure (a) is that for the 
potential energy like that of gravity, for which the graph of U is that in Figure 2(a), puck Shares ee ee 
there will be bound or unbound motion, depending on the energy level. This is Pn WS cea 
illustrated in Figure 2(b), 


yo 


U = $h(r = ra) 


(a) (b) 
Figure 1 


For gravity, the potential 
energy is 
v= aM 


= 
as in Equation (19) of 
Section 2. 

(a) {b) 


Figure 2 


The puck connected to a spring is reconsidered in the light of this analysis. Since 
friction is negligible the force is central, and the effective potential energy looks like 
that of Figure 1(b). From this graph, inner and outer bounding radii (radial turning 
points) are predicted for the motion. Viewing the puck’s motion from an overhead 
camera verifies this prediction. 


In summary, isotropic central forces conserve both angular momentum and total 
mechanical energy. Without solving the equations of motion, the turning values for r 
can be found from knowledge of the potential energy function U(r) and total 
mechanical energy B. 
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Two further physical examples are considered. The first is gravity, which is the main 
topic of this unit. The gravitational potential energy U and its corresponding effective 
potential energy U*" are the functions shown in Figures 2(a) and 2(b) respectively. So 
in this case there may be two turning points, as for elliptical orbits, or only one turning 
point, as for unbound hyperbolic orbits. 


The second example is the mutual repulsion between two electrostatic charges of equal 
sign, either both positive or both negative, Here the potential energy U has the same 
general shape as the gravitational case (Figure 2(a)), but with opposite sign. The 
functions U and U°" for the electrostatic force are shown in Figures 3(a) and 3(b) 
Lelow. It is clear that there are no bound states in this case. The potential U for both 
examples can be written in the form U = ji/r, where j: < 0 for gravity and j: > 0 for 
electrostatic repulsion. 


(v> 0) 


of ’ 


(a) (b) 
Figure'3 


The physicist Ernest Rutherford carried out an experiment in the first decade of the 
twentieth century to demonstrate that atoms consist of relatively massive, extremely 
small nuclei bearing a positive electrostatic charge, surrounded by light electrons 
carrying a compensating negative charge, The experiment consisted of firing the lighter 
nuclei of helium atoms, with some energy F, at other heavier atoms and measuring the 
(perforce) unbound ‘scattering’ orbits of the helium nuclei, Rutherford achieved good 
agreement with his experimental data by applying the classical mechanical theory of 
orbits. Surprisingly, the modern theory of quantum mechanics gives similar theoretical 
predictions, for subtle reasons. 


Summary of Section 4 


1. Ifascalar function @(r) depends upon position only through the radial distance 
|r| =r = (x? + y* + 27)!/?, so that o(r) = d(r), then the gradient of @ can be 
written as 


dd 
rad ¢ = —e,, 
grad¢ = 7 er 
where e, = r/r is the unit vector in the outward radial direction. 
2. Hence if F is an isotropic central force, that is, 
F(r) = f(r)er 


for some function f, then F is conservative, with a potential energy function U(r) 
such that 


The motion of a particle in this force field is planar, and has constant total 
mechanical energy 


EB = }m(i* +16?) + U(r). 
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5 End of unit exercises 


Exercise 1 


For a gravitational orbit, it was shown in Example 1 of Section 3 that the total mechanical 
energy is 


where m is the mass of the orbiting particle, M is the mass of the central body, G is the 
gravitational constant, and e is the eccentricity of the orbit. 


(i) For a bound orbit, show that 


mMG 
2a 
where a is the semi-major axis of the elliptical orbit. 


B=- 


(ii) For any orbit, show that 


p= mG” _ 1), 


2rmin 


where ryyin is the distance of closest approach of the orbiting particle to the origin. 


Exercise 2 


It is proposed to put a geostationary satellite into a circular orbit above the Earth's equator, 
such that the period of the satellite's orbit matches exactly the period of the Earth's rotation. 
Show that the distance d of the satellite from the Earth’s centre is given by 

@ookt 

4?" 

where R is the Earth’s radius and 7’ is the period of its rotation. Taking g = 9.81 ms~? and 
R=6.4 x 10° m, calculate the corresponding height of the satellite above the surface of the 
Barth. 


Exercise 3 


A particle of mass m is launched with speed v at right angles to the line which joins it to a 
spherically symmetric body of much larger mass M. At the instant of launch, the distance 
between m and the centre of M is R. You showed in Exercise 8 of Section 3 that, for the 
subsequent motion, the angular momentum has magnitude 


L. = moR, 

the total mechanical energy is 
$ 2_ GmM 
B= }m?- =", 


and the motion is bound if 
R 
(i) For bound motion, show that the semi-major axis of the orbit is 


aya eee Nh 
o= (3-4) 


(ii) For bound motion, find an expression for the eccentricity ¢ in terms of v, R, M and G. 
(iii) Show that the condition for circular motion is 


2_ MG 
aoe. 


[Solutions on page 52] 


This situation was considered 
previously in Unit 28 
Subsection 4.2. 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) Integration of dU'/de = —k gives U(r) = —kw + C. 
Here C is an arbitrary constant, which may be chosen equal 
to zero for simplicity, leaving U(x) = —ke. 
(ii) Integration of dU/da = ki gives U(x) = thu? + C, 
where C’ is an arbitrary constant. Choosing C’ = 0 produces 
U(a) = th2?. 
(iii) When U(x) = k/.r, the force is 

F(x) =—dU/dx =k/x* (x > 0). 


2. Differentiation of 
a(t) = 29 coswt + ™ sinwt 
gives . 
#(t) = —wa sinwt + v0 coswt. 
(Note here that (0) = ro and (0) = vo, as claimed.) 
Squaring each of these equations gives 


‘ovo 


wx? =a} cos? wt + coswt sinwt + 


and 


wiry sin? wt — Qwiry vp coswt sinwt + vg cos” wt. 
The total mechanical energy is 

Bt) ne Monat whet a CO 

E(t) = gmié° + pha” + C 

= dma? + dma? + C. 

On substituting for x? and #? from above and simplifying, 
this becomes 

B(t) = }mw?x} + dmv} + C = B(0). 
Hence, by direct substitution, it has been shown that F is 
constant, 


At various times the particle reaches its turning points, 
located at a7 say. Since the speed vanishes at these points, 
the energy there can be expressed as 


B= that + C = }mwrey + C. 
Hence, since E is constant throughout the motion, we have 


dmw*ah + C = kmw?xG + dmvg + C, 
from which the locations of the turning points are deduced 
to be 


op = +h + 03/w?. 


3. Since the particle is released from rest, its initial velocity 
(0) is zero, so that its total mechanical energy is then 
E = }mi?(0) + U(a) = Ula). 


The energy is constant throughout. At x = /3 it can be 
written as 


E = 4mw?(3) + U(). 
where u(r) = at. Equating the two expressions for E gives 
Emv?(3) + U(3) = Ula), 


so the speed at w = 3 is 


0(3) = 2 wa) —U()).- 

Similarly, the speed at x = 7 is 
2 

¥(9) = \/ = (Ula) - 0). 


For « > 6 the potential energy function U(:r) is zero, so the 
speed is then constant, with value 


v= 


From Figure 6 of Section 1, or by considering these formulas 
for the speed and remembering that U() < 0,U(y) > 0, we 
have 


v(B) > v(6) > v(y). 


4. The particle is projected from a point to the right of 
x =, where the potential energy vanishes (see Figure 6 of 
Section 1). Since its initial velocity is (0), the energy is 
B= }mi*(0). 
from which we have 
2E 


m 
‘The negative sign here arises because the particle is 
projected to the left. After projection it slows down to a 
minimum speed at x = and then accelerates to an overall 
maximum speed at « = 3. After reaching this point it slows 
to a halt at «=a, Thereafter it retraces its path, with a 
maximum speed at x = @ and minimum speed at «x = ¥. 
After reaching x = 6 it moves out to infinity with velocity 
equal to its initial speed, that is, 

EB 


m 


(0 


t= 


5. Using the definition of the gradient, the vectOr equation 

— grad U, 

with F = —mgj, is equivalent to the three scalar equations 
I OU au 


By inspection, these equations have the solikion 
U=mgy+C, 

where C is any constant. Thus the uniform gravitational field 

is conservative, with the potential energy function just found. 


6. We seek 
F = —gradU = 
where 


U=h/r = ka? + y? +27), 
Now we have 

Ow 

ar 
Similarly, it can be shown that 

au . au 

BU = tyr? ond 


Lk (ant + y? + 27)-9/? x 2 = —kar 


=-ker 
so that 
k F k 
B= S (ai + uj + 2k) = a 


7. Tt is given that 
mt =—gradU +G and G+F=0. (1) 
Differentiating E = }mé + U and applying the chain rule 
leads, as in the text, to the equation 
E=-+(m¥ + gradU). 
Using each of Equations (1) in turn, this becomes 
E=7-G=0. 
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8. As explained in the text prior to the exercise, the total 
mechanical energy E is constant, where 


E=1mu? + mgy+C. 
Initially, the particle’s speed is zero, since it is released from 
rest, and its height is y = h, so that BE = mgh +C. When it 
reaches the origin we have y = 0 and u = uo, so the energy 
can also be written as E = }mug + C. Hence 

dmus + C = mgh+C, 
giving the particle’s speed at the origin as 

uo = /2gh. 
(Note that this outcome is independent, of the slope of the 
plane.) 


9. Since i,j form an orthogonal pair of unit vectors, we have 
i-j=0 and i-i=j u 
‘Then Equations (9) of Section 1 give 
e, + €¢ = (cosdi + sin@j) - (—sin0i + cos@j) 
= —cos@sin i+ i+ (cos* @— sin? @)i-j 
+ sin @ cos Oj -j 
=0. 
Similarly, we obtain 
e, +e, = (cos i + sin8j) - (cos Oi + sin Oj) 
= cos* # +sin?@=1 


and 
eo ey = (sind i + cos@j)+(—sindi+ cos j) 
= sin" 0+ cos’ @ = 1, 


10. Since the particle is constrained to move at a fixed radius 

r= R, we have 7 0. Then, in plane polar coordinates, 
t= Riley and # = —Rée, + Riey. 

(These are the results which you obtained in Exercise 7 of 

Unit 28 Section 4, and which were derived again in 

Subsection 4.2 of that unit.) Using Equations (9) of 

Section 1 to substitute for e,.e in terms of i,j gives 


# = —Résin#i + ROcos Aj, 
# = —R (6? cosd + dsind) i+ R (Gcosd — 6" sind) j, 


11. The polar angle @ is fixed here at the value 0. Hence 
6=6=0. Then, in plane polar coordinates, 
tafe, and #=7e,. 


In terms of the Cartesian unit vectors i and j, these are 
F=f cos Ooi + Fsin Oj, 
cos Ao i + sin A j. 


é 


12. We have 
r=al?+bt, #=2at+b, 
and @=ct*, 6=2ct, 6=2, 
so that the velocity and acceleration are given by 
F = (2at + be, + (at? + bt)(2et)ey 
and # = [2a — (at” + bt)(2ct)"Je, 
+ [2(2at + b)(2ct) + (at? + bt)(2c)Jep. 


2a, 


13. The velocity of the particle is 
i= He, +réey, 
and since e,,eg are orthogonal unit vectors (with e, - ey = 0 
and e, +e, = e+ ep = 1) the total mechanical energy is 
E=imi-t +U 
= 4m (¥* +776?) +U. 


14. (i) The force acting on the particle is F = —Te,. Since 
isa linear function of t, we have # = 0. Hence Newton's 
second law, in the form of Equations (14) of Section 1, gives 
j2 im md / a4 
~mr6?=—T and ™S (r36) 
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These equations are equivalent to 


mré =T, (2) 
d pag 
a (r76) =o. (3) 


(ii) From Equation (3), we know that r°4 is constant. 
Hence we have 


r*(t)(t) = r?(0)6(0). 
This is used to substitute for #(t) in Equation (2), giving 
1?(0)6(0)\* _ mr*(0)62(0) __mr*(0)67(0) 
Pe) ( oo ) > te O ea : ae 
as required. 
(iii) The string will break when 7’ = 500N. From the answer 
to part (ii), this occurs when 


$062 
mr(0)670) _ so 
an 


or when 
r= (mr (0)6(0)/500)'”*. 


With m = 0.1kg, r(0) = 1m and 6(0) = 1rads™*, this gives 
r = (0.1/500)'/* ~ 0.058 m. The time at which the string 
breaks is then 

t=(r(0) = 1r)/n = 60(1 — r) ~ 56.55. 
(Note, from the form of the function for T, that for any finite 


breaking strength the string will break before the particle 
reaches the origin.) 


15. (i) The force is —A(r — Io)e, so that, from 
Equations (14) of Section 1, the equations of motion are 


m (i — 16?) = —k(r — lo), 


d Fy 
ES (v8) 
(ii) From the last equation, mr*d must be constant. (‘This 
quantity is the angular momentum of the particle. 
Conservation of angular momentum occurred also in the 
previous exercise, and will be of great importance in the 
derivation of planetary orbits.) 


m 


(ili) If U = 3h (r — lo), where » = (a* + y?)'/?, then we have 


as tex 2(r = to) 
= h(n — lo) * 3(a? + y?)7/? x Qe 
=k(r—l)ar™*. 
Similarly, 
a =h(r=lo)yr. 
so that (using Equation (8) of Section 1) 
au, ov, 
—grad U = 
= —k(r = lo)r "(ai + yi) 
= =k(r = lo) tr 
= —k(r — Ip)e,. 


Since this expression for — grad U is equal to the force F 
given in the statement of the exercise, we have shown that L/ 
is a potential energy function for F (and hence established 
that F is a conservative force). 
(iv) Suppose that the speed of the particle is denoted by v. 
From part (iii), the total mechanical energy 

E = mv? +U = dmv? + k(n = lo)? 
is constant. When r = 2fo, the speed is given as v = u, so 
that 

B= 4mu? + ih. 
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When r = lo the potential energy is zero, so that the speed v 
is then given by E = }mv?. Hence we have 


Pees ee 
pmu? = pm 


or v=4/u?+ 


m 


for the speed of the puck when the spring has its natural 
length. 


Solutions to the exercises in Section 2 


1. When a = b, the equation of the ellipse can be written as 
X?+Y¥? =a. 
This is the equation of a circle with radius a. 


2. (i) The area A is twice the area of the upper half of 
the ellipse, namely 


A=2 ‘g Y¥(X)dX. 


But from the definition of the ellipse (Equation (1) of 
Section 2), the upper half of the ellipse (where Y > 0) has 


equation Y = by/1— (X/a)?. Thus the area is 
A= » [ V1 (X/a)? dX. 


(ii) Let X 
u=0, and dX 


acosu, When X = -a, u =m, when X =a, 
~asin udu, so that the integral becomes 


A= 20h | sin? udu = aud f sin? udu 
™ Jo 
» at f (1 — cos 2u) du = nab. 
Jo 


3. A circle corresponds to a = b, for which q = 0. Hence, 
for a cirele, the foci coincide at the centre. 

4. (i) When b =a the eccentricity is e = 0. 

(ii) When b = $a, we have « = } V3 ~ 0.866. 

(iii) When 6 = 4a, we have e = 1 V8 ~ 0.943. 


The larger the semi-major axis a is compared to the 
semi-minor axis b, the more ‘squashed’ is the shape of the 
ellipse. ‘This corresponds to larger values of the eccentricity ¢. 


5. From Equation (9) of Section 2, we have 
A=tR'6. 

But also V = RO, so that 
A= $RV. 


Since both R and V are constants, A is constant. 


6. Refer to Figure 4 of Section 2. During the interval 

t; St < ta, the planet is closer to the origin than it is during 
the interval t; < ¢ < ty. So to sweep through the same area 
(as it must do by Kepler's second law), it moves faster 
during the first interval. 


7. IfT is the planet's period and a is its semi-major axis, 
then Kepler's third law states that 


=ka*, 

where k is a constant. If T\r,am are respectively the period 
and semi-major axis for Mars, and Tr, ae are the 
corresponding parameters for the Earth, then we have 


sy = (@y or Tu =Tp (ay. 


For the values given, 
Tm = 365.256 x (1.524)*/? ~ 687.187 days. 
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8. In terms of units, Equation (10) of Section 2 for the 

force of gravity gives 

G x (mass)? 
(length)? 

where ‘=" here means ‘has the same units as’. It follows that 

_ (length)? x force 


force = 


i (mass)? 
But from Newton's second law, F = ma, we have 
force = Mass * length 
(time)? 
so that 
aes (length)? mass x length __(length)® 


(mass)? (time)? ~~ mass x (time)?* 
Thus, in SI units, G has the units m° kg~! s~?. 


9. The given vector is of the form a/la|, where a= rj — r,. 
Any vector of this form is a unit vector, since its magnitude 
is a|/la| =1. 


10, Equation (10) of Section 2 gives 
= Gmumy (9) 

pane 
as the force on m, due to mj. Hence, swapping the indices i 
and j, the force on my, due to m is 

Gmyjm, 

Fy = rar 
where we have used the fact that |r, ~r,|=|ry —ry|. Thus 
F,, = —F);, and Newton's third law is satisfied by the 
gravitational force. 


7 


11. Equation (15) of Section 2 gives the gravitational 
potential energy function resulting from a system of particles 
m,ma,.+..mw, acting upon another particle of mass m 
located at r. Here there is only a single particle in the 
system (N = 1), with mass m, = M and position vector 
r = 0. Hence the potential energy function is 
GmM GmM 
U(r) = - so, 
(r) ry = 
The corresponding gravitational force is 
F(r) = —grad U(r), 
From Equation (13) of Section 2, with r; =r) = 0, this gives 
1 imMr 
F(r) = GmM grad — = 
Ir] {rls 
which could also have been obtained directly from 
Equation (12) of Section 2. Putting e, = r/|r| for the unit 
vector in the radially outward direction, this can be written 
as 


F(x) = - SM 


12. By definition, a spherically symmetric mass distribution 


has a mass density function p which depends on r only 
through its magnitude 


frl=r= Ye? +y? +22. 
Of the functions listed, only (c), (d) and (f) meet this 
criterion, 


13. From Equation (19) of Section 2, the gravitational 
potential energy of the particle of mass m is 


U(r) = 2M, 


This is the same potential energy as was found in Exercise 11 
for gravitational attraction by a single particle of mass M at 
the origin. Hence the gravitational force on the particle of 
mass m, is, as before, 
imM 
B(r) = —gradU(r) = — 2. 
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14. Suppose that Body 1 has mass M and that Body 2 has 
mass m. Using the given hint, the motion of the centre of 
mass C2 of Body 2 (which is also its geometric centre, by 
spherical symmetry) is the same as that of a particle of mass 
m acted upon by the same total force. However, we know 
from the result of Exercise 13 that the gravitational force 
experienced by such a particle due to Body 1 is directed from 
Cz to the centre C; of Body 1, with magnitude GmM/R?. 
This is therefore the total gravitational force exerted by 
Body 1 upon Body 2. 

By a similar argument, or by application of Newton's third 
law, Body 2 exerts upon Body 1 a force of the same 
magnitude but in the opposite direction, that is, from Cy to 
Ci 


15. (i) The magnitude of the gravitational force on a 
particle of mass m at the Earth’s surface (where r = R) is 
mg (from the simple model) and mMG/R? (from Newton's 
universal law of gravitation). Since these must be equal, we 
have g = MG/R?. 

(ii) From the given values, we obtain 

9.81 x (6.4 x 10°)? 


= 24 
Rarerosie = Oe 


16. At height h above the surface, the gravitational force on 
a particle of mass m has magnitude mMG/(R +h)", so that 
the gravitational acceleration at this height is 
MG 

(REA) 

At ground level, we have 
MG 

9= er 
from Exercise 15(i). The gravitational acceleration at height 
his therefore 

OR wrens, 

(R+hP (+ h/RP 
With g = 9.81ms~?, 10m and R = 6.4 x 10° m, this 
gives the value 9.78ms~* for the gravitational acceleration 
at a height of 10km, 


Solutions to the exercises in Section 3 


1, Forces (a), (b), (c) and (f) are central forces, since these 
are directed along r or along —r, depending on the sign of 
the constant K (and, in the case of (c), on the sign of x), 


2. From the result of Exercise 13 in Section 2, the 
gravitational force acting on the particle of mass m is 


F(r) = - SM, 


giving g(r) = -—GmM/r. 
3. (i) By the definition of the cross product of two 
vectors (Unit 14 Subsection 3.5), 
L=rx mt 
is perpendicular to r whenever L # 0. 


(ii) Combining Newton’s second law and the definition of a 
central force (Equation (1) of Section 3), we have 


mi = F(r) = g(r)r, 
for some scalar function g(r). Then we obtain 


er : 
L= F(r x mi) 


W 


PX met+r x mE 

=0+rx o(r)r =0, 
where we have twice used the fact that a x a= for any 
vector a. Hence L is a constant vector. 


(iii) Referring once more to the definition of the cross 
product, if L = 0 and r 4 0,r ¥ O, then ¢ must be in the 
direction of r or of —r. 


MST204 30 Solutions 


4, (i) The angular momentum of the particle is 

L=rx mr, 
but since the motion is in the («, y)-plane, we may apply the 
equations 


r=re, and +=fe,+re9 


(Equations (8) and (12) of Section 1). Since e, x ey =k, we 
obtain 


L=re, x m(re, + rey) 
= mrie, x e, +mr*be, x eo = mrOk. 
Hence L = L.k, where L. = mr*0. 
(ii) Since m is constant, L. = 0 is equivalent to 


d (ag) 
a (r?6) =0, 
which is Equation (3) of Section 3. 


5. (i) Isotropic central forces have the form 

F(r) = f(r)e,. Of the forces listed in Exercise 1, the 
following are isotropic central forces; 

(a) since F(r) = Krer; 

(b) since F(r) = (K/r)er; 

(f) since F(r) = Key. 


(ii) For the gravitational force we have, from Exercise 2 (or 
from Exercise 13 in Section 2), 


80 this is an isotropic central force with f(r) = —GmM/r*. 


6. ‘The effective potential energy function is 
12 GmM 
Qmr? or 
Differentiation of this equation gives 
eS GmM 
a coe 
Now dU*"/dr = 0 at r =r", which leads to the solution 


- ii 


mM 
The corresponding effective potential energy is 


GmM __m'M2G? 


ue = 


vet 


ln = U(r”) 


7. Ife =Usi, then we have # =0 and + =r’ at all times 
This orbit is therefore a circle of radius r*. 


8. (i) The angular momentum is L. = mr?6. At the 
instant of launch we have r = R and, since there is then no 
radial component of velocity, the initial speed is vu = RA. 
Hence L. = mvR. 


(ii) The total mechanical energy # is the sum of the kinetic 
and potential energies, and is constant. Initially, the kinetic 
energy is }mnv? and the potential energy is —GmM/R. 
Hence, for all subsequent times, 


B= Yr? — GM, 
(iii) The motion is bound provided that E < 0. So for 
bound motion, we must have 

2MG 


2 GmM F 
pmo? === <0, or v*< Cre 


9. (i) Since v? = MG/(2R) the total mechanical energy 

is, from Exercise 8(ii), 

2 GmM __3GmM 
0 haa. 

Since £ is negative, the motion is bound. 


=i 
E = bmw 
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(ii) In general, bound motion has two turning points 
(although for the special case of motion in a circle the radius, 
r, is constant). From the given hint, the equation for turning 


points is B= U*"(r), where 
jetty - 22 _ GmM 
OS a a 


From Exercise 8(i), we have L. = muR, so the equation for 
Tmin and rox is 


(moR)* _ GmM _ jp, _ _3GmM 
2mr? 28 4R * 


Putting « = R/r, the equation becomes 
mvs? GmMs __3GmM 
Die Varn ao 

‘The given condition v? = WG/(2R) simplifies this to 
132 i 0. 


4s?-s=-2,. or s*-4s+3 
. which correspond to 


; 
The solutions are s = 1 and s = 3, 
rinax = Rand rin = $R- 


(We must assume here that the spherically symmetric body 


at the origin has a radius less than 4. Otherwise the 
orbiting particle would collide with it!) 


10. (i) The equation of orbits is 
£ = 1+ecos(0~00). (1) 


The right-hand side of this equation takes values between 

nd 1 +e (inclusive) as @ varies. The minimum value of 
¥ corresponds to the maximum value of the right-hand side, 
which is 1+. Hence we have 


nin = tte’ 

(ii) Ife <1 then 1 ~e > 0, so that the right-hand side of 
Equation (1) is always positive, Its minimum value 1 — 
corresponds to the maximum value of r, giving 


Tmax = 


(iii) Ife > 1 then 1~e <0 and the treatment of part (ii) 
cannot apply, since r is always non-negative. The condition 
r > 0 is equivalent to 

1+ ecos(8 — A) > 0. 
which is satisfied by the range of angles 

A — 0, <8 < + Oa. 
where 9, = arceos(—e~'). For this range of values of # there 
is no upper bound for the values of r, so the motion is 
unbound in this case. The corresponding orbits (for ¢ > 1) 
are hyperbolas, as illustrated in the diagram below. 


7’ 


Figure 1 


(iv) Orbits are bound for e < 1 and unbound for e > 1. 
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11. When ¢ = 0 we have (from Equations (18) and (16) of 
Section 3) 
mo 
~ MG" 
This is a circle of radius J, for which rmax = min =. From 
Example 1 of Section 3, the total mechanical energy is 
mMG __m°M?G? 

Ss ane 
(The radius of the circle agrees with the value of r” found in 
Exercise 6, as expected from Exercise 7. Since # = 0, the 
corresponding energy is E = U*"(r*), and the value of E just 
obtained agrees with that for U°"(r*) found in Exercise 6,) 


r= 


B= 


12. (i) From Equation (13) of Section 1, the acceleration is 
#=(#— rb Je, + (216 + rb)ep. 


Here r =d and @ are constants, so that * 6 =0, Hence 


# = -dé*e,, 
showing that the acceleration is in the inward radial 
direction, with magnitude |i| = d0?. 
(ii) From Equations (8) and (12) of Section 1, we have 
ie, + rbey = rbeg, 
so that the angular momentum is 

L=rx mi =mr*de, x eg = mak. 
Thus the magnitude of the angular momentum is L. = 


(iii) The gravitational force on the satellite is (from. 
Exercise 13 of Section 2) 


r=re, and r=* 


md". 


F= one, = Sn e,, 
86 that, from part (i), mit = F gives 
ai? = MG. (2) 
Using the result of part (ii) to eliminate @, this becomes 
ibe 
meade or = TMG" 


in agreement with the value of the orbifal radius { found in 
Exercise 11. 
(iv) The gravitational potential energy is -GmM/d. The 
kinetic energy is (from Equation (2)) 

Lm(dd)? = eat 
so that the sum of the kinetic and potential energies is 

pa Gm 

ene) 

With d = 1, this is the same as the expression for E found in 
Exercise 11. 


13. (i) If the Earth is assumed to be spherically 
symmetric, and d is the distance between the Moon’s centre 
of mass and the Earth's centre, then by Newton's universal 
law of gravitation, the magnitude of the gravitational force 
exerted on the Moon is [F| = GmM/d?. 


(ii) We have || = dé*. The angular velocity @ is constant, 
and therefore equals 27/7’, where T is the orbital period. It 
follows that 
je = 2d 
= 
(iii) The result g = MG/R? is shown as in the solution to 
Exercise 15(i) of Section 2. 


52 


(iv) From Newton's second law m¥ = F applied to the 
gravitational motion of the Moon, and from parts (i) 
and (ii), we have 
4r?7dm _ GmM 
72. ad 
which is equivalent to 


a_ (4) 5s 
4 - (a) * 


From part (iii), this is the same as 


as required, 
(v) Using d= 60R, g = 9.8ms~*, T = 28 x 24 x 3600s and 
R=5.3 x 10° m gives 
Antd’ _ 4x? x 60° x 5.3 x 10° 
gh@T? — 9.8 x 28? x 247 x 36007 
This is the ratio of the two sides of the equation in part (iy), 
which should be close to 1 if the prediction is accurate. (If 
R= 6.4 « 10° m is used instead of 5.3 x 10° m, with the other 
values unchanged, then the value of this ratio becomes 0,95.) 


= 0.8, 


14. (i) Ina time T the whole area of the ellipse is swept 
out. Since A is constant, we then have 


rab = : 
2m 
Squaring both sides of this equation, and putting 
Ey 
# =al= 
oO 7eMG 
we obtain 


ya _ (4) 5 
z = (a rs 


as required, 


(ii) For the circular orbit of Exercise 13, the orbital radius 
is a= b=d, Use of MG = gR? (from Exercise 13(iii)) then 
gives the result of Exercise 13(iv) once more. 


Solutions to the exercises in Section 5 


1, (i) If the orbit is bound, then we have 
> MMG) 2 _ 
BS iceaigl tam 
together with 
v 


t=— and e= 
a 


eaae es (4) and (5) of Section 2). Hence 
2] = —p/q? = 
_mMG 


—l/a, giving 
£ 


2a 

(ii) For any orbit, Exercise 10(i) of Section 3 gives the 
distance of closest approach as 

l 
I+e 
Hence we obtain 
mMG 
——(e-1). 
2rmin ‘ ) 
(As a check. both of these results agree with that of 
Exercise 12(iv) of Section 3 when the orbit is circular, with 
radius d =a =rmin and eccentricity ¢ = 0.) 
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2. The result 
_ ght? 
f= 4n? 


can be obtained just as in parts (i)-(iv) of the solution to 
Exercise 13 in Section 3, though here d is the radius of the 
satellite's orbit rather than that of the Moon's. On putting 
g =98lms~*, R=6.4 x 10° m and T = 1d: 24 x 36005, 
this gives 


= a (9.81 x (6.4 x 10°)? x (24 x 3600)*) 


or d~ 4.24 x 10m ~6.6R. 
‘The corresponding value for the height above the Earth’s 
surface is 
d—R~3.6 x 10" m. 
(The use of geostationary orbits for communications 


satellites was first suggested by the science fiction author 
Arthur ©, Clarke.) 


3. (i) For bound motion, we have both 


2 GmM 
B=im?-S" 
and (from Exercise 1(i) of this section) 
mMG 
baat 


On equating these expressions and solving for a, we obtain 


Pig) GA Ges 
«= (3-aa) ; 


(ii) From Equations (4) and (5) of Section 2, we have 


sfictnyict 


But also, it is known from Equation (16) of Section 3 that 
Be 

m?MG" 

muR is given in the question, we have 


where the magnitude has been taken to ensure that e > 0. 
(Note that the condition e < 1 must correspond to the given 
condition v? <2MG/R for bound motion.) 

Alternatively, this result could have been obtained from the 
formula 


Feel) =i) 


PT] 
(iii) For circular motion, we have ¢ = 0 and 
pe MG 
=> 


